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PREFACE 


This report presents the results of one phase of research carried out 
at the Jet Propulsion Laboratory, California Institute of Technology, for the 
NASA Ames Research Center by agreement with the National Aeronautics 
and Space Administration (Contract NAS 7-100). W. Kendall and A. Ishimaru 
are consultants to the Jet Propulsion Laboratory. 
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Abstract 


The purpose of this study is to predict the effects of turbulence in the 
atmosphere of Venus on Pioneer Venus radio. Estimates of turbulence effects 
depend a great deal on our knowledge of the turbulence. Although very little 
is known about the turbulence in the Venusian atmosphere, valuable information 
can be inferred from radio measurements of the planet's atmosphere made by the 
previous space probes Venera and Mariner 5. 

A careful investigation based on a theoretical and experimental study of 
the power spectrum of the Mariner 5 amplitude fluctuations is carried out and 
the results contribute considerably to our scientific knowledge of turbulence 
in the atmosphere of Venus. Fully developed turbulence is seen to exist pre- 
dominantly in the altitude range of 41 - 49 km. This result is consistent 
with the high wind shear and wind velocities observed by Venera 4 for altitudes 

higher than 40 km. The outer scale size of turbulence is on the order of 

-7 -1/3 

100 m, the structure constant for the dayside atmosphere 3.9 x 10 m , and 

-7 -1/3 

that for the nightside atmosphere 2.9 x 10 m 

Estimates for Pioneer Venus indicate that the variance of log-amplitude 
fluctuations is less than 0.056, the variance of phase fluctuations less than 
0.88 rad 2 and the bandwidth of both amplitude and phase fluctuations no more 
than a few Hz. 
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I. Introduction 


Since the atmosphere of Venus is very dense, even slight perturbations 
in the atmosphere caused by winds, thermal currents, etc., could give rise to 
substantial fluctuations in the refractive index. As a result, a radio wave 
propagating through the Venusian atmosphere would be expected to suffer 
fluctuations and undergo spectral broadening. These turbulence effects play 
a key role in the design of the communications link with the Pioneer Venus 
entry probes and the purpose of this study is to estimate them. 

The quantities that we are interested in predicting for the Pioneer Venus 
link consist of the probability density functions, the means, the variances, 
and the power spectra of the amplitude and phase fluctuations. The communica- 
tions engineer can in turn use these results to estimate the tracking ability 
and degradation of the phase coherent system. To make predictions for the 
turbulence effects we must know the nature and the characteristics of the 
atmospheric turbulence. No direct measurements for the Venusian atmosphere 
have been made. However, some pertinent information can be inferred from the 
radio link of spacecraft that probed the atmosphere of Venus. These of course 
include Mariner 5 and the Venera spacecraft. Our predictions for Pioneer Venus, 
therefore, rely heavily on radio measurements made by previous space probes. 

This report is organized in the following manner. In Section II we present 
a very brief summary of the theory of wave propagation through turbulent medium. 
To include a complete presentation is impossible within the scope of this report. 
Nevertheless, our objective is to point out the salient features and to provide 
the reader with a source of tutorial references. In Section III the basic equa- 
tions are presented and predictions of the effects of Venusian atmospheric 
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turbulence on the Pioneer Venus radio link are given. These estimates are 
based on the Venera radio data. In Section IV we present the results of pro- 
cessing the Mariner 5 occultation data. The turbulence effects that appear to 
be present are then compared with the developed theory as well as the Venera 
results. An additional crude estimate of radio spectral broadening can be 
obtained by analyzing the received spectra from the Venus radar experiment. 
This is discussed in Section V. Finally, conclusions and future work are sum- 
marized in Section VI. 
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II. Wave Propagation Through Turbulent Medium 


In recent years there has been an increased interest in the problem of 
fluctuations in an electromagnetic wave as a result of scattering from random 
variations in the refractive index. In the earth's atmosphere these effects 
are often negligible at S-band, but become dominant in the millimeter and 
optical bands. Consequently, much of the recent work has dealt with the 
higher frequencies. In the case of Venus, because the lower atmosphere is so 
much denser than that of earth, turbulence is expected to affect frequencies 
as low as S-band. The available theory for wave propagation through turbulent 
medium is applicable to both the microwave and optical bands. 

Extensive investigations of turbulence effects on wave propagation have 
been carried out by Tatarski (Refs. 1-2). Although very mathematical, the 
study of wave propagation in a turbulent medium consists of two basic steps. 

The first is the solution of Maxwell's equations to obtain the fluctuations in 
the electromagnetic wave in terms of the refractive index fluctuations. The 
method most widely used is Rytov's method which takes into account diffraction 
and multiple scattering effects. Rytov's method is valid when the fluctuations 
of the electromagnetic wave are small as in the case of microwave propagation 
through the atmosphere of Venus. The second step in the study of wave propaga- 
tion in turbulent medium involves the characterization of the refractive index 
fluctuations. The power spectrum of the refractive index fluctuations is most 
convenient because its form in the atmosphere is the same as that of the wind 
velocity fluctuation spectrum. The wind velocity spectrum has received much 
attention from the atmospheric physicist and there is both theoretical and 
experimental evidence to indicate that its form follows that of the well-known 
Kolmogorov spectrum for fully-developed turbulence (Refs. 1-5, 8, 9). 
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Several excellent review articles of atmospheric turbulence and its 
effect on wave propagation have been published. Two general articles by 
Wheelon in 1959 (Refs. 6, 7) give a review of the state-of-art at that time. 

More recently, Strohbehn (Ref. 8) reviewed Rytov’s method and its range of 
validity, Kolmogorov’s spectrum and the probability density functions of the 
amplitude and phase fluctuations, while Lawrence and Strohbehn (Ref. 9) 
summarized the variances and temporal power spectra for the amplitude and 
phase fluctuations. Although the emphasis of the last two references is on 
optical frequencies, the results are equally applicable to microwave frequencies 
and the validity of this application is discussed in Ref. 10. The state-of-art 
of wave propagation in turbulent medium has also been reviewed in a recent 
Russian paper (Ref. 11) which is particularly useful because it contains a very 
extensive list of references. Several survey papers cover the effects of 
atmospheric turbulence on communications. Brookner (Refs. 12, 13) reviews the 
turbulence effects in terms of a communication channel model, while Kerr, et. al 
(Ref. 14) discuss the consideration of atmospheric turbulence in the design of 
a communication system. 
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III. Estimates for Pioneer Venus Based on Venera Radio Data 


Let us consider a turbulent medium characterized by the refractive index 
n(r) = 1 + n x (r) 

where n-^(r) is the fluctuation and assumed small and r is position. We write 
the propagating electric field in the absence of turbulence 

u 0 <?) - A o( ?) 

where A Q is the amplitude and <}> 0 the phase. When the turbulence effects are 
included the electric field becomes 

i(j> (r ) 

u(r) = A(r) e 

The turbulence effects are generally described in terms of the quantities 

A 

log-amplitude fluctuation x = ^n — and phase fluctuation S = <j> - if> . 

Aq 

We note that (3) can be alternately written as 

u(r) = u G (r) e^ r ^ 

where <p = x + IS. For weak fluctuations. Maxwell's equations can be solved 
using Rytov's method to yield the first iterated solution (Refs. 1, 2, 8) 



where G is the free space Green's function, k is the free space wave 
number and dV' is the differential scattering volume element. In the 
case of Pioneer Venus, we will assume that the radio frequency is 2297 MHz 
so that k = 48.11 rad/m. It is clear that (5) describes the fluctuations 
in the electromagnetic wave in terms of the refractive index fluctuations. 
According to (5) and assuming that the correlation distance of n-^(r) 
is small compared to the dimensions of V' , x and S are both composed 


( 1 ) 


( 2 ) 


(3) 


(4) 


(5) 
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of the sum of a large number of independent components, and therefore have a 
Gaussian distribution by virtue of the Central Limit Theorem. This means that 
A has a log-normal distribution. There have been alternative arguments in 
favor of the Rayleigh or Rice-Nakagami distribution for the amplitude A. For 
weak fluctuations there is little difference between the Rayleigh and log-normal 
distributions. Moreover, experimental evidence seems to support the log-normal 
distribution. A good summary of the distribution of A can be found in Refs. 8, 

9 and 11. 


The variance of the log-amplitude fluctuations is a measure of the 


strength of signal amplitude fluctuations. For weak fluctuations (a^ << 1) : 

A 


<(x - <X>) 2 > 


A - A, 


■o ) 


where the brackets < > denote an average or expectation and <y> - 0. If the 
conservation of power is invoked, we can use the identity (Ref. 15) 


<e ag > 


a < g > +| a 2 <(g 
e 


<g>) 2 > 


(g is Gaussian) 


to show that (4) yields 


<A> 
or <x> 



As mentioned in Section II, the refractive index fluctuations are 
characterized by the power spectrum which we denote $ n (K), where K is the 
spatial wave number. For well-developed turbulence, $ n (K) follows the 
well-known Kolmogorov spectrum, which can be written (Refs. 1, 2, 8, 9) 
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$ n 00 


(9) 


= 0.033 c 2 (r) K -11/3 


where the structure constant c n characterizes the intensity of turbulence 
and may vary as a function of position. The turbulence present in the 
atmosphere of Venus appears to be well-developed (Refs. 16, 17) and we will 
use Kolmogorov's spectrum of the form proposed by von Karman (Refs. 1, 2, 8) 
to describe it, i.e.. 


*n<K) 


0.033 c 2 (r) K + 


1 - 11/6 


L 


where L 0 is the outer scale of turbulence. The reason for using (10) instead 
of (9) is that (10) includes the effects of blob sizes larger than L Q . 
Additional comments on this point will be made later. 

The spatial correlation functions for X and S can be formed in terms 
of the spectrum (10) using (5), and the respective variances obtained. In 
a previous study (Ref. 18) it was shown that when the path length in the 
turbulent atmosphere was shorter than kL 2 , as in the case of Pioneer Venus, 
then the effect of outer scale size on a 2 is negligible. For the Pioneer 
Venus configuration shown in Fig. 1, we therefore have (Refs. 1, 2, 17) 


(z = z Q ) 


0.565 


c n(z) 


5/6 

z dz 


( 10 ) 


(ID 


where a 2 (z = z Q ) is the variance of the log-amplitude fluctuations of the 
radio signal received from the entry probe when it is located at altitude z c . 
As can be seen, the spatial variation of c n plays an important role in the 
determination of the effects of turbulence. Meaningful prediction of signal 
fading for Pioneer Venus due to atmospheric turbulence cannot be made without 
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some reliable estimates of the structure constant and its variation with 
altitude. Although no direct measurements of the structure constant have 
ever been performed, estimates can be inferred from the radio data trans- 
mitted through the Venusian atmosphere by space probes. In this section we 
will discuss the Venera data only. The Mariner 5 results will be covered in 
Section IV. 

Veneras 4 through 8 were slow-descent probes with Veneras 4, 5 and 6 
reaching approximate altitudes of 26, 24 and 19 kms respectively, and 
Veneras 7 and 8 successfully landing on the surface of Venus. The signals 
received from the Venera spacecraft were processed and their standard 
deviations computed over various time segments as the spacecraft descended 
through the atmosphere (Refs. 19-21). The standard deviations due to 
instrument error and interplanetary plasma scintillation were estimated 
and subtracted from the processed results, and the remaining part a x was 
attributed to atmospheric turbulence. Fig. 2 is reproduced from Ref. 21 
and summarizes the results obtained for Veneras 5-7. To estimate the 
structure constant, Yakovlev, et. al. (Ref. 21) assumed that c n was constant 
in the altitude range 0-55 km. In other words, they assumed that the 
turbulence was homogeneous for altitudes lower than 55 km and non-existent 
for higher altitudes. Thus, (11) can be easily integrated to yield 

c = 0.308 c£ k 7/6 L 11/6 (12) 

where L is the penetration distance of the probe into the Venusian atmosphere. 

-7 -1/3 

The dashed line in Fig. 2 represents (12) where a value of 1.4 x 10 cm for 
°n waS use< ^' Although it is obvious that there is a large scatter of data, it 
should be pointed out that much of this could be due to the variation of c n 
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with the various missions. A more meaningful check would be established 
if just the Venera 7 results were used. Unfortunately, Yakovlev, et. al., 
did not label their data. Although there does appear to be less scatter 
when (12) is compared with the separate sets of data obtained from Veneras 4 
through 6 (Refs. 19, 20), the estimates derived from the Venera radio data 
are at best order-of-magnitude estimates only. 

Prior to Venera 7, de Wolf (Refs. 22, 23) extrapolated the Venera 4 
results to the surface of Venus and pointed out that serious fading could be 
expected there. This extrapolation assumed that the structure constant 
increased exponentially with decrease in altitude and was reasonable since 
this type of variation has been observed in the earth's lower atmosphere 
(Ref. 24). On the other hand, the homogeneous model used by Yakovlev, et. al . 
is also plausible in view of the fact that although the lower atmosphere is 
dense, it is relatively calm as evidenced by the atmospheric wind velocity 
profiles (Ref. 16). At any rate, the difference between extrapolating the 
structure constant inferred from Venera 4 exponentially or as a constant, is 
significant. In the case of exponential extrapolation, severe fading could 
occur. In the absence of radio data from the surface of Venus, other than 
that of Venera 7, we will assume that the turbulence is homogeneous. It must 
be emphasized that our predictions are only as valid as the estimates of 
turbulence inferred by Yakovlev, et. al., from the Venera 7 radio data. 

According to Yakovlev, et. al. (Ref. 21), when Venera 7 was on the 
surface of Venus a was 0.073. For Venera 7, X = 0.3 m and letting L = 55 km, 

A 

! 

(12) yields 


"n 


1.012 x 10 -6 nf 1/3 


JPL Technical Memorandum 33-644 


9 



It might be interesting to note that (13) is about an order of magnitude 
higher than the value for c n deep in the earth's atmosphere (Refs. 24, 25, 28). 
Using the value of c n in (13), we obtain from (12) the following estimates for 
Pioneer Venus on the surface of Venus: 


o x = 0.1189 

o2 = 0.014 


( 14 ) 


and the corresponding <A> can be computed using (8) . 

In the case of phase, = a |. If we assume (10), we can derive the 
following results for a | in homogeneous turbulence using the general results 
given in Ref. 18: 


o| = 0.782 k 2 L 0 5/3 c 2 L - a2 (15) 

where o^ is (12). As is expected, the phase fluctuations, unlike the 
amplitude fluctuations, receive contributions from blob sizes larger than 
Lo (Refs. 8, 18). Although L 0 is unknown, some estimate can be derived 
from the Mariner 5 data as will be shown in Section IV. L appears to be 
approximately 100 m. If L q is 100 m, we let L = 55 km and from (12) and (15) 
we obtain for the probe on the surface of Venus 


= 0.4535 rad 

aj = 0.2057 rad 2 


(16) 


It should be pointed out that according to Ref. 20, Venera 4 observed stronger 
turbulence than Veneras 5—7. In fact, the value of c n is approximately twice 
that inferred from Veneras 5-7. If we use this worst case estimate, we have 
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0.2378 


al = 0.056 

A 

Oq = 0.9373 rad 

= 0.8786 rad 2 


(17) 


Until now we have discussed only the variances of the log-amplitude 
and phase fluctuations. It is important in the design of a coherent 
communications link to know how these fluctuations are distributed in 
frequency. We shall therefore consider the temporal frequency spectra of 
the log-amplitude fluctuations W^(f) and the phase fluctuations Wg(f). If 
the time correlation functions of x and S are B x (t) and B s (t), respectively, 
their Fourier Transforms are the corresponding frequency spectra. Thus 


oo 

W x (f) = 4 | B x (t) e" i2lTft dt 

OO 

W s (f) = 4 j B s (t) e" i2irft dt 

— OO 

To determine the variations in time and consequently the frequency spectra 
of the log-amplitude and phase fluctuations, Taylor's frozen- in hypothesis 

is generally assumed (Refs. 1, 2, 9). Under this assumption, the atmospheric 

inhomogeneities are frozen-in and the time variations of the received signal 
are due to movement of the inhomogeneities across the propagation path by the 
wind. Essentially, only the component of wind velocity perpendicular to the 
propagation path need be considered and we denote it v . The frozen-in model 
is valid as long as the lives of the individual eddies of the turbulent 
atmosphere are longer than the time it takes for an eddy to cross the propa- 
gation path. Tatarski (Ref. 2) has shown that this condition is met when 
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L << L 2 /A. For L = 100 m, L 2 /A = 77 km. The frozen- in hypothesis is, 
o o o 

therefore, valid even when the Pioneer Venus probe is on the surface of Venus. 
Another point worth mentioning is the fact that Tatarski (Refs. 2, 9) has also 
analyzed the case where the motion of the eddies is completely random. His 
results indicate that the spectra for completely random motion and those for 
constant motion (frozen turbulence) are similar. 

We assume that the Pioneer Venus probe in Fig. 1 is transmitting a 
spherical wave. Then, using Rytov's method and the frozen-in hypothesis, 
the received spectra at a distance L' away are (Ref. 25): 


where 


w x (f) 

w s (f) 


L’ 





h Y (n,K) 

= 8TT 2 k 2 

K dK 

dn 

A 


. 


h s (n,K) 


J J 

2uf 0 


V K > 


1 

/k 2 vj^ - ( 2 tt f ) 2 


h x (n,K) 

h s (n,K) 


1 - 

1 + 


n , v k 2 
cos - (L — n ) - 
l n _ 

cos -jr (L - ri) 


(19) 


( 20 ) 

( 21 ) 


Letting L' = L + L" (see Fig. 1) and noting that L" >> L, L 1 in (19) becomes L, 
and (20) and (21) reduce to 

h x (n,K) = 1 - cos 

hg(n,K) = 1 + cos 


( 22 ) 

( 23 ) 
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It should be noted that (22) and (23) along with resultant (19) are the same as 
the frequency spectra for a plane wave received at a distance L away. Again 
we will use $ n (K) as given by (10) and will assume the turbulence to be 
homogeneous. In addition, we will assume that the distribution of vj^ 
throughout the atmosphere is constant even though the Venera probes indicate 
that it is far from it (Ref. 16). Still, the wind velocity distribution 
should not significantly affect the spectra in terms of their bandwidth, a 
parameter which we are primarily interested in estimating. In the future 
we hope to evaluate (19), (22) and (23) with the aid of a computer taking 
into account the wind velocity profiles observed by the Venera space probes. 


When L 0 

= the following 

asymptotic expressions 

can 

be obtained (Refs. 1, 


2, 9) 

w x (f) = 

f-K) 

°X 

0.44 ~ 
± 

where 

f J. 

. i li 

2ir /L 

(24) 


w s (f) = 

2W x (f) 




(25) 


f-K) 

f-x» 






w x (f) = 

f-x» 

32.26 W x (0) 

ff r 8/3 

KJ 

= 

Oy f 'I —8/ 3 
1 13 * 

jl rj 

(26) 


w s (f) = 

K X< £ > 




(26) 


f-K» 

f-X» 





It must 

be emphasized that 

(24) - (26) are asymptotic values only and the 



fact that Wg is written in terms of and vice versa does not mean that 
and Wg are dependent on each other. 
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The computer was used to calculate the intermediate values as well as the 

spectra corresponding to finite values of L Q . These results are summarized in 

Figs. 3 and 4 and are normalized with respect to W (0) for L 0 = i.e., (24). 

i /iT ^ 

The various values of 3 = — j— plotted correspond to L Q = 50, 100 and °° m 
for L = 55 km. The dependence of the spectra on L 0 indicates the consequences 
of using von Karman's representation of the Kolmogorov spectrum (10) instead 
of (9). For the amplitude fluctuations, the variation of the spectra with L c 
is small. The reason of course is that the contribution of scattering from 
blob sizes larger than L Q to the amplitude fluctuations is very small. On the 
contrary, phase fluctuations do depend on scattering from blob sizes larger 
than L 0 and the difference between von Karman's representation (finite L Q ) and 
the Kolmogorov spectrum (infinite L 0 ) is consequently most pronounced at 
frequencies lower than f r Another reason for using von Karman's representation 
is demonstrated in Fig. 4 where it can be seen that the Kolmogorov spectrum 
(L 0 = 00 ) leads to a singularity in Wg(f) at the origin. 

Many features of the spectra shown in Figs. 3 and 4 have been observed 
by experiment at various microwave and laser frequencies in the earth's 
atmosphere (Refs. 26-29). The most consistent is the -8/3 slope at the higher 
fluctuation frequencies, which is a consequence and therefore verification of 
the Kolmogorov spectrum. The breakpoint in W (f) which occurs at approximately 

A 

f “ f has also been used to derive estimates of c n and v (Refs. 25, 28). 

As previously mentioned, an important parameter in the frequency spectra 
for the communications engineer is the bandwidth. Denoting the frequency 
corresponding to the half power point for L 0 = «=, f^g, we have 
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We note that neither the shape nor the half power point of the log-amplitude 
spectrum is dependent on c n or a^. On the other hand, f^ dB is proportional 
to v, . It is also clear from Fig. 3 that the dependence of bandwidth on L 0 
is very small. This of course means that the choice of either (9) or (10) 
as the spectrum has little effect on the bandwidth. The spectra given in 
Figs. 3 and 4 and (27) should be valid for L >> /XS, Q , where £ 0 is the inner 
scale of turbulence. In the earth's atmosphere, l Q is of the order of 1 cm. 

If we assume a similar value for Venus, the spectra and (27) should be valid 
for distances as short as 0.036 m. An interesting result of (27) is that the 
spectrum becomes narrower as the entry probe penetrates deeper into the 
Venusian atmosphere. Actually, this is not surprising in view of the fact 
that at a distance L, the predominant scale size is the Fresnel-zone size Al" 
(Refs. 9, 18). Thus, for longer L, the predominant scale size increases, and 
the contribution to the lower frequencies in the frequency spectrum correspond- 
ingly increases since v^ remains constant. To summarize, as L increases the 
frequency spectrum becomes narrower even though the variance increases. 

We have obtained estimates for Pioneer Venus based on the values of 
structure constant and inferred from Venera 4 and these are shown in 
Table 1. Venera 4 is the worst case analysis since the peak horizontal wind 
velocity and structure constant inferred from Venera 4 were twice those 
inferred from Venera 7. Additional calculations of f 3dB are summarized in 
Fig. 5. 

The results discussed so far assume that the sub-earth direction corre- 
sponds to the zenith of Venus (Fig. 1). If the sub-earth direction makes an 
angle 0 with the zenith of Venus, then L in (12), (15) and (27) must be replaced 
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by L sec 0, while in (27) must be replaced by cos 0. These results are, 

of course, not valid when 0 is large (0 = 90°). We see that while cr 2 an d a 2 
increase with 0, the spectra become narrower. Thus as far as the bandwidth is 
concerned, 0=0 represents the worst case. 

In conclusion, let us once again emphasize that our predictions for 
Pioneer Venus rely heavily on our knowledge of the dependence of the structure 
constant on altitude and the measurement and calculation of a 2 from the Venera 
radio data. Let us further point out that a 2 has a very small value and is 

A 

consequently subject to substantial experimental error. More important, however, 
is the fact that we have little assurance that the fluctuations are due to 
turbulence and not noise or some other cause such as absorption layers in the 
atmosphere or spacecraft motion. Computation of the frequency spectrum in addi- 
tion to a 2 and its ' comparison with the spectrum shown in Fig. 3 would go a long 
X 

way toward providing supporting evidence for atmospheric turbulence. Although 
some form of spectral analysis was performed on the Venera data (Ref. 30), its 
purpose was not to compare the spectrum with theory. There are two reasons for 
raising this point. First of all, we will demonstrate the importance and use- 
fulness of the spectrum in our analysis of the Mariner 5 data in the following 
section. Secondly, it should be apparent from our discussion that with proper 
analysis, the Pioneer Venus radio link could be utilized to provide valuable 
information on turbulence in the atmosphere of Venus. 
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IV. Observations of Turbulence In the Atmosphere of Venus by Mariner 5 


It should be clear from the previous section that in order to make 
engineering estimates for Pioneer Venus, we must have information on the 
turbulence characteristics of the atmosphere of Venus and how it interacts 
with wave propagation. On the other hand, the only way we are going to gain 
this knowledge is through interpretation of radio signals received from pre- 
vious missions to Venus. The engineering and scientific aspects of atmospheric 
turbulence are therefore closely related, and consideration of one cannot be 
made without the other. It is obvious that in spite of the Venera missions, 
our knowledge about turbulence in the atmosphere of Venus is severely limited. 
The confidence in our predictions for Pioneer Venus would be enhanced if 
similar results could be obtained based on the radio data of Mariner 5. 

However, since Mariner 5 was a fly-by mission, new techniques for analysis 
and data interpretation will have to be developed. 

Let us first review what has been done with the Mariner 5 data to date. 

In his study, Gurvich (Ref. 17) computed the variance of log-amplitude fluc- 
tuations 0^ for Mariner 5 during occultation using the published results of 

A 

Kliore, et. al. (Ref. 31). The structure constant inferred from 0 ^ was found 
to be at least two orders of magnitude lower than that inferred from Venera 4. 
De Wolf (Refs. 22, 23) also investigated Mariner 5, and concluded that the 
observed scintillations were much weaker than those that would have been 
predicted from the observed Venera 4 scintillations because the Mariner 5 
fluctuations had been smoothed by data processing. Golitsyn and Gurvich 
(Ref. 32) examined the feasibility of extracting the turbulence-induced 
fluctuations from the noisy Mariner 5 exit occultation radio signal, and 
concluded that the signal-to-noise ratio was sufficiently high in the S-band 
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link. In a recent paper, Woo and Ishimaru (Ref. 18) showed that the analysis 
employed by Gurvich (Ref. 17), although adequate for entry missions, was 
inadequate for fly-by missions because it neglected the effects of turbulence 
inhomogeneity and finite extent in the direction transverse to the direction 
of propagation. A new analysis for was developed and it was shown that 

the structure constant inferred from Mariner 5 was within an order of magnitude 
of that inferred from Venera 4. 

In a totally different interpretation, Rasool (Ref. 33) suggests that the 
rapid changes in the Mariner 5 amplitude data may be due to absorption by cloud 
layers containing substantial amounts of Hg and its compounds. Thus, even the 
presence of turbulence effects in the Mariner 5 data is questionable. 

It is clear that a careful and thorough study of the amplitude and phase 
fluctuations of Mariner 5 is highly desirable. This study should be based on 
the power spectra of the amplitude and phase fluctuations instead of just the 
variances. Two approaches should be used: a theoretical analysis of Mariner 5 

and spectral analysis of the Mariner 5 occultation data. It should be pointed 
out that the power spectra and phase fluctuations have never been studied in 
previous work. We have undertaken such a study of Mariner 5 and its purpose is 
(1) to examine and ascertain the effects of data processing, (2) to distinguish 
turbulence- induced fluctuations from those due to noise, (3) to obtain supporting 
evidence of atmospheric turbulence in the atmosphere of Venus, (4) to derive the 
theoretical spectra for the turbulence- induced fluctuations of Mariner 5, (5) to 
infer the turbulence characteristics of the Venusian atmosphere by comparing the 
processed and theoretical spectra, and (6) to compare the inferred results of 
Mariner 5 with those of Venera 4 - 7 in order to enhance our confidence in our 
estimates for Pioneer Venus. The results obtained so far and described in this 
report are those based on the amplitude fluctuations. 
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4.1 Mariner 5 Processed Results 


Although the Mariner 5 occultation data were processed in previous work 
(Refs. 31, 34 and 35), the primary interest was the average atmosphere so that 
substantial smoothing was generally introduced into the data and little attention 
was paid to the very delicate measurement of variance. Since we are interested 
in the fluctuations, we have had to reprocess the Mariner 5 data, a somewhat 
difficult task in view of the age and volume of the data. 

A complete description of data processing would make this section prohib- 
tively long. We have therefore described most of the details in Appendix A, 
which is organized in the following manner. In A.l signal reception and pre- 
processing are discussed. The pre-processing includes digitization, frequency 
shifting, filtering and decimation of the original Mariner 5 analog tapes. In 
the present study, the starting point is the digitized tapes. The phase- locked- 
loop processing is described in A. 2 and consists of a discussion of amplitude 
and phase extraction, errors and tracking. Finally in A. 3 we present an account 
of the amplitude spectrum calculation. 

As mentioned earlier, the results obtained to date deal only with the 
amplitude fluctuations. Although they are preliminary, they are nevertheless 
significant. The signal level as a function of time during entrance and exit 
occultation is shown in Figs. 6 and 7, respectively. These results are similar 
to those previously published (Refs. 31 and 32). We have also indicated the 
altitudes corresponding to the lowest penetration of the propagation path. The 
information contained in Ref. 35 was used for the case of entrance occultation 
while that in Ref. 32 for the case of exit occultation. We expect to obtain 
more accurate data on the exit case in the future. We have divided Figs. 6 and 
7 into three regions: A, B and C. The reason for this is that the fluctuations 
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in Figs. 6 and 7 suggest noise in regions A and C and turbulence in region B 
which corresponds to the approximate altitude range of 41-49 km. Of course, 
verification of these results will emerge from the spectral study. As discussed 
in Appendix A, we fit quadratic time functions to Figs. 6 and 7 assuming that they 
represent the average atmospheric attenuation as well as the changing antenna 
pointing and earth-spacecraft distance effects. The quadratic time function is 
then subtracted from the data in Figs. 6 and 7 and the results spectrum analyzed. 

The power spectra of the amplitude fluctuations for the various regions 
are shown in Figs. 8-12. At the writing of this report only the results for 
region A of the entrance case had not been completed. The spectra are plotted 
for frequencies up to 1/2 ttt^q£ since this is the maximum frequency at which 
amplitude fluctuations may be observed (see Appendix A) . Noise levels as well 
as 90% confidence intervals were calculated using the formulas developed in 
Appendix A and are shown on the various spectra. The 90% confidence interval 
is indicated by error bars. Errors of the spectral estimates are expected to 
lie within the error bars in 90% of all cases. For smoothing, we used running 
point averages consisting of 11, 21, 31, 41 and 51 points (m is the number of 
points) . 

Let us first examine the results obtained for the exit case. It is 
apparent that the spectra of regions A and C do indeed resemble noise (flat 
spectrum), and that in contrast the spectrum for region B shows structure above 
the noise level at the lower frequencies. The calculated noise levels seem to 
agree with the apparent noise levels in the spectra. The spectra for the 
entrance case are similar to those for the exit case. In region B the apparent 
noise level appears to be twice as high as the calculated noise level. The 
important point, however, is that again the spectrum in region B has structure. 
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A careful examination of the results for region B reveals the difficulties 
involved in spectral analyzing the Mariner 5 data. On one hand, because the 
spacecraft is moving very fast, the record length available for spectral analysis 
is very short. As a result, there is a severe limitation on the low frequency 
end of the spectrum. On the other hand, the signal-to-noise ratio is so low that 
the spectrum very quickly reaches the noise level at the high frequency end of 
the spectrum. As a result, frequency averaging would make the already very 
narrow frequency range of the spectrum even narrower. Fortunately, despite these 
limitations, the results are good enough to allow the distinction and separation 
of noise and structure. 


We have utilized the computer to calculate the probability density function 
A 

of x = In — t the log-amplitude fluctuations. The result for Region B in the 
A o 

entrance case is shown in Fig. 13. It can be seen that the probability density 
function for X resembles a gaussian function, so that the probability density 
function for A would be log-normal as expected if the fluctuations were turbulence- 
induced (see Section III). It should be pointed out that for small fluctuations, 
the difference between log-normal and gaussian functions would be small. The 
result in Fig. 13 is consequently not too meaningful. 

Although the structure in the spectra of region B is very likely due to 
turbulence-induced fluctuations, we definitely need additional supporting evidence. 
We hope to provide part of this with the results of our theoretical study which 
will be discussed in the next section. For the moment, we will assume that the 
structure in the spectra in region B represents turbulence. The processed results, 
therefore, indicate that turbulence exists predominantly in the altitude range of 
41-49 km in the atmosphere of Venus. It is interesting to point out that strong 
wind shear, a condition necessary for the existence of turbulence, was indeed 


observed at altitudes higher than 40 km by Venera 4 (Ref. 16) and more recently 
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by Venera 8 (Ref. 37). Kerzhanovich, et. al. (Ref. 16) have recently discussed 
the presence of turbulence in the altitude range 45- 53 km in view of the wind 
velocity observations of Venera 4. The processed results of Mariner 5 are 
therefore significant because they furnish for the first time information on 
the variation of turbulence in the atmosphere of Venus as a function of altitude. 


4.2 Mariner 5 Theoretical Results 

The results obtained in Section 4.1 indicate that turbulence in the 
atmosphere exists predominantly in the altitude range of 41-49 km. We will 
therefore assume that the structure constant c n has the following spatial 
variation 


(x, 


c’ y c* z c ; 



(28) 


where the coordinate system is depicted in Fig. 14 and the peak structure 
constant c no occurs at the altitude of 45 km. The extent of the turbulence 
in the x c - z c plane is characterized by a, while that in the y c direction 
by b. These two characteristic distances are approximately related by the 
following relationship (see Fig. 14): 

a = 6095 tan0 (29) 


-1 

cos 


e 


6095 

6095 + b 


where a and b are in kilometers and the radius of the planet Venus is assumed 
to be 6050 km. Since the turbulence exists in the altitude range of 41-49 km, 
b is expected to lie in the range 4-6 km. The corresponding range for a is 
220 - 270 km. It should be pointed out that (28) assumes essentially no turbulence 
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for altitudes less than 35 km. Of course no information was provided by 
Mariner 5 since it did not probe below this altitude on account of super- 
ref ractivity . We will discuss this point further later in the report. 

Since the outer scale size L c is expected to be somewhere around 100 m, 
the distance L (see Fig. 14) which is of the order of 10 4 km substantially 

o 

exceeds kL Q . The theory for wave propagation in homogeneous turbulence 
utilized in Section III for analyzing entry missions such as Venera and 
Pioneer Venus is therefore no longer adequate (Refs . 18 and 36) . To account 
for the effects of turbulence inhomogeneity and finite extent in the direction 
transverse to the direction of propagation, the planetary atmosphere is 
represented by a localized homogeneous and isotropic turbulence with smoothly- 
varying mean characteristics. Although a general study of the variances of 
log-amplitude and phase fluctuations for wave propagation in localized 
turbulence has been conducted (Ref. 36), no study of the corresponding frequency 
spectra has ever been made. Furthermore, the results described in Refs. 18 and 
36 do not include the effects of the moving spacecraft. The mathematical 
details of the present theoretical study are rather involved and consequently 
described in Appendix B. The formulation of the general problem of wave 
propagation in localized smoothly-varying turbulence is given in B.l. The 
variances are derived in B.2 and the frequency spectra in B.3 and B.4. 

For the sake of comparison we will normalize the frequency spectra of 
the log-amplitude and phase fluctuations, W^(f) and Wg(f), respectively, to 
their values at f = 0. W^(0) and Wg(0) can be obtained from (91)- (95) in 

Appendix B: 
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(30) 


W x .(0) = K'(W’-W’) 

W s (0) = K'(W[+W’) (31) 

where 

l 0 8/3 

K = 0.033 c 2 Q 4 it 3 — k 2 a exp 



Some typical normalized spectra of the log-amplitude fluctuations are 

plotted in Figs. 15-17. The most prevalent result is that the major con- 

1 v 

tribution to the spectra occurs at a frequency approximately equal to YU TT * 
In other words, the predominant scale size is L 0 . This is reasonable since 
the predominant scale size is the Fresnel— zone size /XiT when L << kL Q . 

When L >> kL 2 , as is the case in Mariner 5, the largest scale size is L Q 
and this then becomes the predominant scale size. 

The spectra in Figs. 15-17 correspond to the range of parameters 
expected for entrance and exit occultation of Mariner 5. Several features 
are apparent. First of all, there is little difference between the spectra 
for entrance and exit occultation. Secondly, the shape of the spectrum 
depends to a large extent on the values of b and L Q . Thirdly, the shape of 
the spectrum is not dependent on velocity v, at least for the range of values 
considered (100- 150 m/sec). In the next section we will compare the 
theoretical spectra with the processed spectra of Mariner 5 and infer the 
turbulence characteristics of the Venusian atmosphere. 
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4.3 Turbulence Characteristics of the Atmosphere of Venus 

We have compared the theoretical spectra in Figs. 15 and 16 with the 
processed spectra of Mariner 5 for region B (Figs. 9 and 11). We selected 
a best fit theoretical curve by adjusting only the ordinate. The best fit 
for both entrance and exit occultation occurs when L Q = 100 m, b = 6 km and 
the corresponding value of a is 270 km. We have superimposed this best fit 
on the spectra with varying amounts of frequency averaging in Figs. 9 and 11 
and the agreement is remarkable. It can be seen that the theoretical spectra 
serve very well to define the shape of the processed spectra. 

The processed spectra indicated the existence of turbulence in the 
approximate altitude range of 41-49 km. The y c dependence of c n is plotted 
in Fig. 18 and it can be seen that the value of 6 km obtained for b in the 
best fit theoretical curve is consistent with the altitude variation obtained 
in Section 4.1. 

If we examine the theoretical spectra plotted in Figs. 15 and 16 we see 

that they depend on v, L, b, a and Lq. We know v and L to a good degree of 

accuracy from the trajectory information and b and a to a reasonably good 

degree from the data in Figs. 6 and 7 and the processed spectra. As a result, 

the process of matching the theoretical and processed spectra essentially 

provides a valuable means for measuring L 0 . In the case of Pioneer Venus, 

the phase fluctuations depend very much on scattering from large eddies of 

turbulence. The outer scale size L Q defines the largest scale size of 

turbulence and consequently is an important parameter for predicting the 

phase fluctuations of Pioneer Venus (see equation 15 in Section II). A 

discussion is given in Ref. 16 pointing out that the expected L Q for the 

Venusian atmosphere would be larger than that of the earth's atmosphere. 

In the earth's atmosphere L 0 generally lies in the range of 10 - 50 m. The 

correlation function of the fluctuations of wind velocity inferred from 
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Venera 4 also indicates that L Q would measure 100- 200 m (Ref. 16). Thus, the 
estimate of L Q derived from Mariner 5 is consistent with these aforementioned 
results . 

2 

Let us infer the peak structure constant c no from o x< For entrance 
occultation was 1.856 x 10 while for exit occultation was 

3.398 x 10 -2 . It comes as no surprise that these results are considerably 
higher than those obtained by Golitsyn and Gurvich (Ref. 32). It is obvious 
that too much smoothing was used in the data studied by Golitsyn and Gurvich, 
so that the contribution from the higher frequencies was lost. This of course 
points up the necessity of a spectral analysis even if one is only interested 
in the variance and supports our concern that the measurement of variance due 
to turbulence in the case of Venera may be inaccurate. 

The appropriate region for the values obtained for L Q , b and a is 
region II (see Fig. B-4 in Appendix B) and we will use (49) in Appendix B 
to infer c nQ . For the nightside atmosphere. 


"no 


2.9 x 10 


r 7 -1/3 


m 


while for the dayside atmosphere. 


'no 


3.9 x 10 


-7 - 1/3 


m 


As would be expected, the dayside atmosphere is more turbulent than the 
nightside atmosphere. However, since the difference in temperature between 
the dayside and nightside atmospheres is small (Ref. 35), the differences 
in c no would also be expected to be small. 

In summary, the results obtained from Mariner 5 consistently support 
the assumption that the fluctuations observed are indeed turbulence-induced. 

It should also be pointed out that since Kolmogorov's spectrum was used in our 


(32) 


(33) 


26 


JPL, Technical Memorandum 33-644 



theoretical spectra, the good agreement with the processed spectra provides 
evidence that the turbulence is indeed fully developed. In our future work 
we hope to study the effect of the shape of the spectrum of the refractive 

index fluctuations on W (f). 

X 


4.4 Comparison of Mariner 5 and Venera Results 

As mentioned in the previous section, the wind velocity profiles measured 
by Venera 4 showed high wind velocities and wind shear above 40 km and low wind 
velocities below. The processed results of Mariner 5 indicate that turbulence 
begins to disappear for altitudes lower than 40 km and are consistent with the 
Venera wind velocity observations. Thus, even though Mariner 5 only probed to 
an altitude of 35 km, it is reasonable to expect that the turbulence would 
continue to be absent or remain weak from 35 km down to the surface of the 
planet. The lack of turbulence in the lower atmosphere could be explained by 
the fact that large amounts of energy that might not be available are required 
to move the very dense atmosphere. 

Our predictions for Pioneer Venus in Section III were based on a model 
of the turbulent atmosphere characterized by c n being constant from 0 to 55 km 
and zero above 55 km. Let us re-examine . these predictions in view of the 
Mariner 5 findings. We will assume the variation of c n shown in Fig. 19; 
c r is constant (c^ = c no ) in the altitude range 41 - 49 km and zero elsewhere. 
Using (11) we obtain the following result: 


0.308 k 7/6 L 2 11/6 


11/61 


1 - 


'no 


(34) 
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When L l = 0, (34) reduces to (12). As in Section III, we will infer c no from 

the measurement of by Venera 7 (o^ = 5.3 x 10 ) . For — 41 km and 

X X 

L2 = 49 km, we obtain 


c no 


2.1 x 


10- 6 B - 1/3 


Comparing (35) with (13) we see that the difference is small. The predictions 
made for Pioneer Venus in Section III therefore remain representative. 

Comparing (35) with (32) and (33), we find that c n inferred from Mariner 5 
is within an order of magnitude of that inferred by Venera 7. The important 
point however is that confidence in our predictions for Pioneer Venus has been 
greatly enhanced as a result of the Mariner 5 study . 


( 35 ) 
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V . Venus Radar Experiment 

An additional crude estimate of radio spectral broadening due to 
turbulence in the Venusian atmosphere can be made obtained by analyzing 
the received spectra from the Venus radar experiment. 

During his radar scan of Venus, Goldstein (Ref. 38) reported that the 
observed spectral broadening between the half power points was 0.07 Hz. We 
can use (27) to obtain a crude estimate of v^ by taking L to be twice the 
extent of the atmosphere, i.e., 110 km. We obtain v^ = 0.9 m/sec which is 
within an order of magnitude of the average v^ observed on Venera 7 (Ref. 16). 


JPL Technical Memorandum 33-644 


29 



VI . Summary and Conclusions 


We have shown that our ability to determine the effects of turbulence 
in the atmosphere of Venus on Pioneer Venus radio depends a great deal on 
our knowledge of the turbulence. Although very little is known about the 
turbulence in the Venusian atmosphere, valuable information can be inferred 
from the radio link of spacecraft that have probed the planet's atmosphere. 

These of course include the Venera and Mariner 5 spacecraft. Our predictions 
for Pioneer Venus, therefore, rely heavily on radio measurements of the 
Venusian atmosphere made by previous space probes. 

Even though Venera 7 reached the surface of Venus, the information 
provided on atmospheric turbulence was very limited. Still, we estimated 
the effects of turbulence on Pioneer Venus based on Venera 7, but pointed 
out the urgent need for additional turbulence information to establish 
confidence in the estimates. 

In an effort to increase our understanding of turbulence in the atmos- 
phere of Venus, we have conducted a thorough and careful investigation of the 
Mariner 5 occultation data. The results described in this report were obtained 
by studying the amplitude fluctuations. Unlike previous work, the present 
study is based on the power spectrum rather than the variance of the amplitude 
fluctuations. The results obtained so far have considerably advanced our scien- 
tific knowledge of turbulence in the atmosphere of Venus. These results were 
made possible by the development of new techniques in two major areas: data 

processing and theory of wave propagation in turbulent medium. The study of 
data processing revealed its effect on the amplitude fluctuations and enabled 
noise to be distinguished from the turbulence-induced fluctuations. On the 
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other hand, the theoretical study of wave propagation in localized smoothly 
varying turbulence led to the theoretical prediction of the variance and power 
spectrum of the Mariner 5 amplitude fluctuations. 

The results obtained with Mariner 5 can be summarized as follows. 

Mariner 5 provides evidence of turbulence in the atmosphere of Venus, and for 

the first time, information on its variation with altitude down to 35 km. The 

remarkable agreement between the processed and theoretical power spectra of the 

amplitude fluctuations furnishes strong evidence that the fluctuations are indeed 

turbulence- induced. Turbulence exists predominantly in the 41 - 49 km altitude 

range, and appears to be fully developed since Kolmogorov's spectrum was used to 

describe the refractive index fluctuations. The variation of turbulence with 

altitude is consistent with the high wind velocities and wind shear observed by 

Venera 4 above 40 km. The outer scale size of turbulence L 0 is on the order of 

100 m. The inferred structure constant for the dayside atmosphere is 

3.9 x 10 -7 m -1 ^ 3 and, as expected, it is slightly lower for the nightside atmos- 
-7 - 1 /3 

phere (2.9 x 10 m ' ). 

It is shown that the structure constants inferred from Mariner 5 are 
within an order of magnitude of that inferred from Venera 7. The estimates 
made for Pioneer Venus radio which were based on Venera 7, therefore, remain 
valid. The significant result of our study of Mariner 5 is that the confidence 
in our estimates for Pioneer Venus has been greatly enhanced. 

Our estimates for Pioneer Venus indicate that the variance of log-amplitude 
fluctuations should be less than 0.056, the variance of phase fluctuations less 
than 0.88 rad 2 , and the bandwidth of both amplitude and phase fluctuations no 
more than a few Hz. 
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In future work we plan on (1) reducing and inferring turbulence from the 
Mariner 5 phase data, (2) studying the dependence of the Pioneer Venus power 
spectra of amplitude and phase fluctuations on wind velocity distribution and 
spectrum of refractive index fluctuations and (3) studying the cross correlation 
between amplitude and phase fluctuations for Pioneer Venus. 
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TABLE 1 



- 50 m/sec and uniform (Venera 4 data) 

2 2 

Note thata^, o q and O^a re independent of 


Table 1. Summary of Log -Amplitude and Phase Fluctuation Results 
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Figure 1. Pioneer Venus Entry Configuration 
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Figure 2. (J-% for Veneras 5-7 as a Function of Altitude. 

Taken from Reference 21. 
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Figure 3. Power Spectra of Log-Amplitude Fluctuations 
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Figure 4. Power Spectra of Phase Fluctuations 
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Figure 7. Signal Level vs Time in UT for Mariner 5 Venus Exit 
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Figure 8c. Region A of Exit Occultation: B = 5 Ha, t _ = 0- 33 secs, 
h = 0. 023 secs, m = 21 
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Figure 8f. Region A of Exit Occultation: B = 5 Ha, t = 0. 33 secs, 
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Figure 9d. Region B of Exit Occultation: B = 5 Hz, r = 0. 033 secs 
h = 0.033 secs, m = 31. AGC 
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Figure 9f. Region B of Exit Occultation: B = 5 Ha, T ~ 0- 033 secs, 
h = 0. 033 secs, m = 51. 
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Figure 10b. Region C of Exit Occultation: B = 5 Hz, t = 0. 033 secs, 
L _ . _ „ -Lj AuU 

h = 0.03secs,m=ll. 
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Figure lOe. Region C of Exit Occultation: B = 5 Hz, T = 0. 033 secs, 
h = 0. 03 secs, m = 41. L AGC 
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Figure 11c. Region B of Entrance Occultation: B = 5 Hz, r = 0. 033 secs, 
h = 0. 03 secs, m = 21. L AGC 
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Figure lid. Region B of Entrance Occultation: B = 5 Hz, t = 0. 033 secs, 
h = 0 03 secs, m = 31. ^ AGG 
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Figure !2a. Region C of Entrance Occultation: B = 5 Ha, T = 0.033 secs, 

L# AGC 

h =0.03 secs, m = 1. 
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Figure 12b. Region C of Entrance Occultation: B T = 5 H*. 7 - ^ = 0. 033 secs, 
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h = 0. 03 secs, m = 1 1 . 
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Figure 12c. Region C of Entrance Occultation: B T =5 Hz, 

Li AVjU 
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Figure 12d. Region C of Entrance Occultation: B - 5 Hz, t = 0-033 secs, 
h = 0.03 secs, m = 31. L AGC 
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Figure 12e. Region C of Entrance Occultation: B =5 Hz, T = 0. 033 secs, 
h = 0. 03 secs, m = 41. AGC 
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Figure 12f. Region C of Entrance Occultation: B =5 Hz, r 

h = 0. 03 secs, m = 51. AGC 


- 0. 033 secs. 
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TURBULENT ATMOSPHERE 



Figure 14. Mariner 5 Occultation Configuration 
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Figure 15. Theoretical Power Spectrum of the Log -Amplitude Fluctuations 

for Entrance Occultation for Various Values of L_, b and a: 

v = 150 m/sec, t =48 secs, L = 8, 400 km. 
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Figure 16. Theoretical Power Spectrum of the Log -Amplitude Fluctuations for 

Exit Occultation for Various Values of L^, b and a: v = 150 m/sec, 

t =48 secs, L = 12, 000 km. 
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Figure 17. Theoretical Power Spectrum of the Log -Amplitude Fluctuations for 

Exit Occultation for Various Values of v : Lq = 100 m, b = 6 km, 

a = 270 km, t =48 secs, L = 12, 000 km. 
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ALTITUDE, km 


Figure 18. Variation of Normalized Structure Constant With Altitude 
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APPENDIX A 


SIGNAL PROCESSING TO EXTRACT TURBULENCE-INDUCED SIGNAL FLUCTUATIONS 

In this Appendix we will describe the processing steps which are required to 
first extract the amplitude and phase variations which are imposed on a signal 
received from a spacecraft occulted by Venus, and to then measure the spectra of 
these variations and separate the resulting spectra into components attributable to 
(1) atmospheric turbulence, and (2) noise at the input to the receiver; 

A. 1 SIGNAL RECEPTION AND PRE-PROCESSING 
A. 1.1 The Received Signal 

The signal transmitted from the spacecraft at time t will be received on 
earth some time later due to delays which can be attributed to several different 
effects. Generally these effects can be classified as Venusian atmospheric 
effects, and trajectory effects. By trajectory effects we mean simply the delay 
component caused by the propagation time (at the speed of light) required for the 
signal to propagate in a straight line from the spacecraft to the earth. For our 
purposes this component can be assumed known. There are two atmospheric effects 
which cause a further signal delay, viz., the lengthening of the propagation path 
due to refraction around Venus, and the slower propagation speed through the 
Venusian atmosphere. Both of these atmospheric effects can be divided into two 
more effects, viz., (1) the delay caused by the "average", or non-turbulent , 
Venusian atmosphere, and (2) deviations from this average caused by turbulence 
and other atmospheric non— homogeneities. It is this last component which it is 
desired to analyze. 


JPL Technical Memorandum 33-644 


A- 1 



Let the delay which the signal received at time t has undergone be denoted 


by d(t). [The signal received at time t was transmitted from the spacecraft at 

time t-d(t)]. Then we can write 

d (t) = d (t) + d (t) + d (t) (1) 

tan 

where 

d^_(t) = delay caused by the spacecraft trajectory, 
d (t) = delay caused by the "average" Venusian atmosphere, 

Si 

d^(t) = delay caused by non-homogeneities in the Venusian atmosphere. 

The signal which is transmitted from the spacecraft for processing on the 
ground is a phase-modulated carrier with the carrier not completely suppressed. 
For estimating the signal delay the (assumed constant) carrier at frequency f 

s c 

will be processed. Therefore, the signal received at time t is simply 

s (t) = A s/2 sin{2ir f [ t— d (t) ]} (2) 

s c 

where 

A = rms received signal amplitude. 

The phase (in radians) of this signal is 

$(t) = 2irf [t-d(t) ] (3) 

sc 

and the frequency (in Hz) is 


f(t) 


^ = f 
2tt sc 


f d'(t) 
sc 


= f 

sc 


- f [d* + d' + 
sc t a 



(4) 


A. 1.2 The Ground Receiver 

For our purposes a simplified model of the ground receiver, as shown in 
Fig. A-l, can be used. In effect the receiver simply adds (unavoidable) noise 
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to the received signal, bandpass filters the sum, multiplies the result by a 

constant-frequency local oscillator (LO) which shifts the signal down to base 

band and also up to an image frequency, and then lowpass filters the result to 

reject the high-frequency image. The result of all this is just to add receiver 

noise to the signal and to shift all frequencies down by the LO frequency f . 

LU 

Therefore, the output of the receiver can be represented as 

r (t) = n(t) + A sin {2 ttI <f - f )t - f d(t)]}. (5) 

sc LO sc 

This signal, along with a reference timing signal, is recorded on an analog tape 
recorder for subsequent processing. 

A. 1.3 Digitization 

The next step in signal processing is digitizing the data for computer 
processing. If a signal is bandlimited to a bandwidth W Hz, then it can be 
represented exactly by its samples spaced no more than 1/2W sec apart. Math- 
ematically, this can be expressed as 



i 


where 


. . sin TTX 
sinc(x) = 

TTX 


r = r(t ) = samples of the received signal, 
i i 

t = time of i^ equispaced sample r., 
i 1 


At = time between samples. 
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to the bandpass tin d lowpass falters in the re cexver, the signal ^-Ct) is 
bandlimited to less than 80kHz, so it could be sampled at any rate greater than 
160,000 samples/sec without losing any signal information. The actual signal is 
sampled at a rate of 320,000 samples/sec. This produces a digitized signal which 
contains all spectral components up to 160 kHz. 

A. 1.4 Frequency Shifting 

The spectrum of the (digitized) recorded signal r(t) is typically as shown 
in Fig. A-2a. It consists of broadband noise plus a signal component at the time- 
varying frequency f (t) — f . The purpose of the filter— and— decimate operation 
is to first frequency shift the signal so that the time variation of the signal 
frequency is eliminated, and then to pass the result through a very narrowband 
filter, and, finally, to reduce the data sampling rate to a value corresponding 
to the new narrow bandwidth. Once this is done, the signal will be represented 
by relatively few samples which can be efficiently further processed in a 
computer. 

The signal component of the received signal r(t) is 

r (t) = A 4~2 sin {2ir[(f - f )t - f d(t)]}. (7) 

s sc LO sc 

This is multiplied by the output of a programmed local oscillator whose frequency 
is f (t) , i.e., it is multiplied by -J~2 sin[2irf (t)tj. The result is 

s ^ 

A sin {2 tt [(f - f + f (t))t - f d(t)]} 

sc LO S sc 

+ A sin{2iT [( f - f__ - f c (t))t - f d(t)]}, (8) 

sc LO S sc 
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i.e., half of the signal power is shifted up in frequency by f (t) and half is 

s 

shifted down by the same amount. In practice the frequency f (t) is taken as a 

s 

spline-curve which compensates for the trajectory and average-atmosphere com- 
ponents of the signal's frequency shift. In other words, f (t) is taken to be 

s 

f (t) = 80 kHz - {f - f - f (d* + d’)} (9) 

£> 1 sc LO sc t a J 

where 

f = carrier frequency transmitted from spacecraft, 

f = receiver local-oscillator frequency, 

d^ = time derivative of trajectory delay component, 

d' = time derivative of average atmosphere delay component. 

Si 

Therefore, the upshifted signal component appears at frequency 

f - f n + f(t) - f d'(t) = 80 kHz - f d'. (10) 

sc LO S sc sc n 

Of course, in order to do this it is necessary for the frequencies f and f 

sc LO 

and the delay-rate components d’ and d' to be known. At this point we can 

t a 

accomplish this by simply defining these quantities to be our estimates of them. 
This will result in any errors in these estimates being included in the delay 
component d^(t) which, up until now, we have attributed to atmospheric non- 
homogeneities. Thus it must be remembered that in later stages of the processing 
these errors in our estimates of the other delays must be recognized and removed. 

This will be done primarily be assuming that the errors in our estimates if d' 

t 

and d' are more slowly varying functions of time than is d', i.e., separation 
a n 

will be done in the spectral domain. 
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The noise component of the received signal r(t) is shifted both up and down 
(half the power in each direction) in frequency by f (t) , just as the signal 
component is. The total upshifted spectrum is therefore as shown in Fig. A-2b, 
and the downshifted spectrum is as shown in Fig. A-2c. The spectrum of the 
multiplier is then the sum of the spectra shown in Figs. A-2b and A-2c. Note, 
however, that in the vicinity of 80 kHz only the upshifted spectrum of Fig. A-2b 
contributes to the sum. Therefore, a narrowband filter used to pass only the 
spectral components near 80 kHz will pass the signal components (half of the 
signal power), and half of the noise power present before multiplication. Thus, 
at no loss in signal- to-noise ratio, we have converted the received signal into 
a signal whose phase is 

2w[(f - f + f (t))t - f d(t)J = 2Tr[8xl0 4 t - f d J + k, (11) 

scLOS sc scni 

where k^ is an arbitrary (and unimportant) constant phase offset. The atmospheric- 

induced delay d will be extracted by further processing of this signal, 
n 

A. 1.5 Filtering and Decimation 

As indicated above, the output of the multiplier (mixer), which performs the 
multiplication of the received signal by the programmed local oscillator, is 
passed through a narrowband filter centered at 80 kHz. If the bandwidth of this 
filter is A , then the spectrum at its output will be similar to that illustrated 
in Fig. A-3. 

It has already been indicated that a signal which is band limited to a 
bandwidth W can be completely represented by only its samples taken every 1/2W 
sec. More exactly, this theorem (the "sampling theorem") states that if a signal 
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is sampled every 1 / 2 W sec, then its spectral components in each of the frequency 
intervals iW to (i+l)W, i=0,l,2,***, are "folded" into the interval 0 to W. 

Of course, if the signal has no spectral components above W, then no higher 
frequency components are folded down, and the original signal is completely 
represented by its samples. However, if higher-frequency components are present, 
then a component at frequency f in the interval iW < f < (i+l)W is "folded" down 
to appear at frequency f t| , where 

f - iW i even 

(i+ 1 ) W - f i odd. ( 12 ) 

Thus, the frequency intervals corresponding to even i are simply frequency 

shifted down to the baseband (0,W) , and the frequency intervals corresponding to 

odd i are shifted down to the baseband, but frequency reversed. This effect is 

used to shift the spectrum shown in Fig. A-3 down to the baseband (0,A), by 

reducing the sampling rate of the signal. 

The original sampling rate of the recorded signal r(t) is 320,000 samples/ 

sec. This is reduced by a factor of N (the signal is decimated by N) by simply 

th 

keeping only every N sample output from the narrowband filter of bandwidth A . 

For this, N is taken to be odd, and the bandwidth A is adjusted to be 
A = 160,000/N. The new sampling period is then N/ 320, 000, so in the sampling 
theorem we have 

Sampling Period - ^ - 55^535 (13) 
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or 


W = 


160,000 

N 


A. 


(14) 


Thus, after decimation all the signal energy will appear in the baseband (0, A). 
Furthermore, all signal and noise components before decimation lie in the fre- 
quency interval 

80,000 - A/2 < f < 80,000 + A/2. (15) 

But A is chosen so that 

80,000 = NA/2. (16) 

Therefore, this signal-frequency interval can be expressed as 


NA/2 - A/2 < f < NA/2 + A/2 


N-l 

2 


A < 



A 


or, since A = W and the decimation factor N is chosen to be odd, 

iW < f < (i+l)W 


(17) 

(18) 


(19) 


where 



Thus, the decimation parameters are chosen so that the signal and noise all lie 
in one frequency interval, so that after decimation they will appear in the base- 
band (0,A) without any frequency "foldover" or aliasing. 

Normally N is chosen so that i is even (a value of N = 101 was used for the 
results contained in this report) , so the signal of bandwidth A centered at 80 
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kHz is simply downshifted in frequency by an amount I (N— 1) / 2 J A = 80 kHz - A/2 

due to the decimation. (The amount of decimation which could be used is limited 

only by how small a bandwidth A will contain the signal frequency component 

f d'.) Note that the frequency of the signal after this operation is 
sc n ' 

A/2 - f d' , and the phase of the signal is 
sc n 

4 = 2ir[(A/2)t - f d ] + k (20) 

sen 1 

The signal and noise amplitudes, and their fluctuations, have, of course, been 
preserved through all of these pre-processing steps. 

A. 2 PHASE-LOCKED-LOOP PROCESSING 

The signal v(t) from the receiving and pre-processing operations is given by 
v(t) = A(t)v / 2 sin{27r[(A/2)t - f d (t) ] + k } + n(t) (21) 

where now we explicitly indicate that the amplitude A(t) is a function of time. 
This time variability is due to four factors: (1) the increasing distance from 

the spacecraft to the receiver, (2) the changing orientation of the spacecraft 
antenna as the spacecraft passes behind the planet, (3) the changing average 
atmospheric attenuation as the received signal passes deeper and deeper through 
the planetary atmosphere, and (4) the variable turbulence-induced scattering of 
the signal as it passes through the planetary atmosphere. It is this last compo- 
nent which is of primary interest here. 

In order to quantitatively determine the effects of planetary turbulence, 
it is necessary to extract (or estimate) the amplitude A(t) and the phase ^^^(0 
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from the signal v(t). This is accomplished by computer processing which extracts 


the phase with a digital phase-locked loop (PLL) and the amplitude with a digital 
coherent amplitude estimator. 

A. 2.1 Phase Extraction 

Here we will not give a complete description of the workings of the phase- 
locked loop, since several good references on this subject (Ref s. A1-A3) are readily 
available, instead, we will simply note that the PLL is a non=lincar phase- 
tracking device which can normally (when it is tracking) be well approximated as 
a device whose input is of the form 


x(t) = A/2 sin{2irf t + 0., (t)} + n(t) 

vco 1 


( 22 ) 


where f ^ is the quiescent frequency of the PLL's VCO, and which extracts the 
phase 0^(t) from the signal. In the process of doing this, the PLL passes the 
phase 0^ and the noise n(t) divided by the amplitude A, through a lowpass filter 
to produce 0^ a noisy estimate of the phase 0^. In other words, the output 
phase estimate 0^(t) is that which would be obtained by passing the input phase 
deviation 0^(t) through a lowpass filter, plus an error which can be thought of 
as coming from passing white noise through the same lowpass filter. 

When using the PLL to estimate the signal delay component due to atmospheric 
non-homogeneities, the VCO quiescent frequency is set at A/2 so that 


0, (t) = -2irf d + k . (23) 

1 sc n 1 
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The phase estimate 0 ^ ( t) is then simply noise plus a term proportional to the 
desired delay component d^ lowpass filtered. 

The transfer function of the equivalent lowpass filter of the linearized 
second-order PLL is (Ref. Al, Eq. 5-19) 



1 + ( 25 /W ) s 

n 

1 + (2?/W ) s + s 2 /W 2 
n n 


(24) 


where 

B = one-sided loop noise bandwidth (Hz), 
Li 

R = Tausworthe's damping parameter, 

W = loop natural radian frequency 
n 

£ = loop damping factor. 

Note that these parameters are related by 


_ 4s 2 + 1 

B L 8? W n ’ 



(25) 


(26) 


W 

n 


4B 

Li 

R+l 


(27) 
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The magnitude and the squared magnitude of the filter transfer function L(j2irf) 
are shown in Figures A-4a and A-4b, respectively. 

A. 2. 2 Phase Errors 

The phase estimate e^Ct) produced by the PLL will differ from the desired 
phase 0^ = 2fff gc d + k , due to two factors, viz., (1) higher-frequency variations 
in 0^(t) will be filtered out by the lowpass filter L(s) (unless the loop band- 
width is made sufficiently large), and (2) the low-frequency components of the 
noise n(t)/A, which in this case has an essentially flat spectrum, will pass 
through the filter L(s). Thus, it is clear that for the best possible phase- 
estimate accuracy, the loop bandwidth used should be as small as possible, con- 
sistent with the requirement that it pass all significant variations in the input 
phase 0^. 

If we let 


4> = Phase error (rad) 


= e x (t) - 0 2 (t) 


and 


S (f) = one-sided spectral density of <j> due to noise, 


(29) 


(30) 
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then for spectrally flat noise, which is typical of the receiver noise present at 
the input to the PLL, with one-sided power spectral density watts/Hz, we get 


N n 2 

s (f) = — r | L(j2TTf) I 

* A 2 


(31) 


2 

where N^/A is the spectral density of the noise at the input to the equivalent 

2 

lowpass filter L(s), and A is the power in the signal into the PLL (see Ref. A2, 

Sec. 2.8). Thus, the spectral density of the noise-induced phase error will have 

the shape shown in Fig. A-4b, and a level which depends on the signal-to-noise 
2 

ratio (SNR) N^/A . 

2 

The mean squared value a. of the phase error is just 

9 



S (f) df 
9 


- r 

‘i 


|L(j2wf) | df 


(32) 


2 

A 

In order for the linearized model of the PLL which we have been assuming to be 

2 

valid, this mean squared phase error must be less than about O.l(rad) . Thus, 
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for all the PLL processing described in this report the loop bandwidth was 
chosen small enough to satisfy 


N B 

2 0 L 

= — o 1 °‘ 1 

4 1 * 2 


( 33 ) 


B 


L ~ 10 N 


0 


SNk. in i-Hz Banawiatn 


10 


(34) 


A. 2. 3 Phase Tracking 

As indicated above, one requirement on the choice of loop bandwidth B is 

Lt 

that it be large enough to pass the significant variations in the input phase 
0^(t). An important aspect of this is the requirement that the frequency rate of 
the signal at the input to the PLL not be too large to follow. This is because 
the second-order PLL tracks a frequency rate with a phase offset of (Ref. A2, 
Table 2.1) 


♦ 


ss 


ttR (R+l) 
8 B 2 R 

Li 


(35) 


where 

R = frequency rate (Hz/sec), 

B = loop bandwidth, 

Li 

R = Tausworthe’s damping parameter. 
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The value of the damping parameter R used was R = 2, 


so <p can be expressed as 
ss 


4> 


ss 


ttR 



9_ 

8 



(36) 


For the PLL to track, this steady-state phase error must be kept below about 
90° = tt/ 2. Therefore, the bandwidth must satisfy 



(37) 


or 


B > i/R 

Lt 


(38) 


which, when combined with the phase-error requirement above, leads to 


/R < B < 

i-i 


10 N 


0 


(39) 


An obvious possibility when tracking signals transmitted from deep space by an 
accelerating spacecraft is that 


10 N 


/ R 


0 


(40) 


due to a low SNR and a high doppler. In this case reliable PLL tracking would 
not be possible. It was precisely to counteract this situation, by removing most 
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of the signal's frequency variations, that the pre-processing described in 
Section A.l was performed. 

A.2.4 Amplitude Extraction 

The estimates of the signal's time-varying amplitudes A(t) which are used to 
detect and characterize the turbulence of the Venusian atmosphere are produced by 
a coherent amplitude-estimating circuit which works in conjunction with the PLL. 
Actually, two such circuits are implemented as part of the PLL. One of these acts 
as an AGC circuit for the PLL, and its parameters are chosen to optimize its per- 
formance in that role. The other is nearly identical, but it is used only to 
produce amplitude estimates for further processing, and its parameters are chosen 
to obtain good performance in that role. 

A block diagram of the coherent amplitude-estimating circuit, as implemented 
in computer software, is shown in Fig. A-5. The input signal is v(t) from the 
pre-processing operations described in Section A.l: 

v(t) = A(t) /2 sin{27r[ (A/2) t - f d (t) ] + k } + n(t) (41) 

sc n 1 

where 

A(t) = time-varying signal rms amplitude, 

A = signal-plus-noise bandwidth, 

f = frequency of transmitted signal, 

sc 

d (t) = delay caused by non-homogeneities in the Venusian atmosphere. 
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This can be simplified notationally to 


v(t) = A(t) /l sin[wt + 0^(t)] + n(t) 


(42) 


by letting 


03 = tcA 

= 2 7T f 

vco 

where 

f = quiescent PLL VCO frequency 
vco 

and, as in Section A. 2.1, by letting 


(43) 


0 (t) = -2irf d + k. . (44) 

1 sc n 1 

To further simplify the notation we will drop the explicit time dependence of 
A(t) and 0^(t), denoting then by simply A and 0^. This results in the form for 
v(t) shown in Fig. A-5. 

As shown in Fig. A-5, the input signal v(t) is first divided by the "gain" 

W , which is derived by the rest of the circuit. The circuit is designed to 
g 

attempt to adjust W so that the ratio A/Wg is held at unity, so that the rms 

g 

signal amplitude into the PLL is constant at unity. In the amplitude-estimating 

circuit used only for producing estimates of A(t), the amplitude-stabilized signal 

v(t)/W is not fed to the PLL or any other circuit. It is used only in the 
g 

amplitude-estimating circuit itself. In the amplitude-estimating circuit used as 
an AGC for the PLL, the signal is fed to the PLL as shown in Fig. A-5. 
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A nearly coherent signal of the form sin (cot + 0 ) is obtained from the PLL 

as described in Section A. 2,1, and mixed with the signal v(t)/W to produce the 

g 

signal which drives the AGC circuit. The phase 0^ is the phase estimate of 0^ 
produced by the PLL, and 0^ - 0^ is the PLL tracking error. The output of the 
mixer can be divided into several components. First, the signal is mixed into 
a low-frequency component (A/W^) cos(0^ - 0^) , and a double-frequency component 
at radian frequency 2oo, which is subsequently rejected. The noise into the mixer 
is of the form 

n i (t) n ? (t) 

— = — sin (cot + 0 ) + cos (cot + 0 ) (45) 

W W 2 w 2 

S g g 

where n^(t) and n^(t) are orthogonal components of the input noise, each of which 
has the same spectral density as the input noise n(t). Therefore, the noise 
component in the mixer output is 


^(t) 


W 

g 


SI 


. 2 
sm 


n (t) 

(u)t + e 2 ) + — 

g 


Si sin (cot + 0^) cos (cot + 0^) 


n x (t) 

Si w 

g 


+ Double-Frequency Components. 


(46) 


The output of the mixer has unity subtracted from it, and is then passed 
through an integrator which effectively rejects the double-frequency components, 
and, along with the gain A , determines the responsiveness of the loop. Finally, 
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the signal is exponentiated to determine the loop gain W , and the estimated rms 

g 

signal amplitude is taken to be W . 

g 

The performance of the AGC loop described here was first analyzed in Ref. A4. 
This analysis was then generalized by Tausworthe in Chapter 7 of Ref. Al. The loop 
analyzed by Tausworthe is shown in Fig. A-6, which is taken from Figure 7-1 of 
Ref. Al. It is shown in Refs. Al and A4 that the loop can be very well approximated 
as a linear loop if the input and output' are identified as the input signal power 
and the estimated signal power in decibels . To this end, following Ref. Al let 


a(t) = 20 log A(t) 

= input signal power in dB 


(47) 


and 

a* (t) = 20 log W 

g 

(48) 

= estimated signal power in dB . 

Then the AGC circuit can be represented as approximately linear with a loop 

i 

transfer function C(s) between a(t) and |c(t) given by (see Ref. Al, eq. 7-6) 


C (s) 


K 


AGC 


[1 + K 


AGC 


Y(s) 

Y(s)] K 

A 


(49) 


where 


K e K 
AGC 

^AGC 20 log e 
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The quantities in these equations can be related to the parameters used in the 
actual processing as shown in Fig. A- 5 as follows. 


Y (s) 


K. 


K 


= 20 , 
= 1, 

= A 


log e 


T 1 n 1 0 T 

AGC AGC 


X AGC = ^" n P ut P arair - St er to the processing software. 


These parameter values give 


K. 


AGC 


AGC 


and 


(50) 


C (s) 


1/s T 


AGC 


0.05 


11 + 1/st agc ! 20 


1 + ST 


AGC 


(51) 


which result in the input signal power (in dB) a(t) being related to the estimated 
* 

signal power a (t) by the block diagram shown in Fig. A-7. Note that the input 
to the equivalent RC filter (represented by K C(s)) consists of three terms: 

(1) the true signal amplitude a(t) in dB, 

(2) a bias of 20 log cos <}> (t) due to phase tracking errors in the PLL, 

(3) a noise component [20 (log e) n, (t)]/W . 

1 g 

Note also that the noise component depends on the value of the gain W . However, 

s 

in most cases of interest the SNR will be high enough for the error in W to be 


A -20 
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at most moderate, so that W 
approximated by 


g 


A, 


and the noise component can be adequately 


20 (log e) n (t) 

7 — . (52) 

A 

The power spectral density of n^ is the same as that of the input noise n(f), 
viz., Nq. Therefore, the spectral density of the noise into the equivalent 
circuit shown in Fig. A-7 is approximately 


(20 log e) N 


0 



(53) 


where the term in brackets is the input reciprocal SNR in a one-Hz bandwidth. 
This flat spectrum will be shaped by the squared magnitude of the equivalent RC 
filter shown in Fig. A-7. The result will be an output error spectrum having 
the shape shown in Fig. A-8. 

A. 2. 5 Amplitude Errors 

There are three basic sources of error in the amplitude estimate a (t) ; one 

for each of the input terms shown in Fig. A-7. These are: 

(1) variations in the true amplitude a(t) are filtered by the RC filter 

with time constant t of Fig. A-7, and if these variations are not 

AGC 

slow compared to t , the estimate a (t) will not follow the rapid 
fluctuation components (which suggests using a small value for t ), 
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( 2 ) 


the noise term proportional to n^t), filtered by the noise bandwidth 

& 

of the filter, is added to the estimate a (t) (which suggests using a 
large value for x ) , 

(3) the phase-error term -20 log cos 4> (t) , which is always positive, is 

•Jg 

filtered and subtracted from a (t) , biasing it toward estimated values 
which are low. 

Here we will treat only the second two of these error sources, and will assume 
that the bandwidth of the equivalent RC filter is large enough (x is small 
enough) to pass the true amplitude fluctuations of interest. 

A. 2. 5.1 Bias Errors 

Since the low-frequency response of the filter in Fig. A-7 is unity, the 
expected value of the estimate a*(t) is simply 

a*(t) = a(t) + 20 log e n (t) + 20 log cos <f> (t) (54) 

x/2 W 1 

g 


Now we can use the fact that, for reasonably small phase tracking errors 

2 

<Kt) = 0^ - 0^, the approximation In cos <j> » -<(» /2 is valid, to write 

2 

log cos <Kt) = - (log e) a / 2 


(55) 


where a is the PLL mean-squared phase error. At high SNR's it is well approx- 

2 2 

imated by (see Section A. 2.2) a = N B /A , where B is the PLL noise bandwidth. 

U L Li 


Also, the average value of the noise n^(t) is zero, so we have 
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( 56 ) 


a*(t) = a(t) - (10 log e) a 
= I(t) -4.34 N q B l /A 2 


The term A /N B is the SNR within the PLL bandwidth, and for our application it 
0 L 

will exceed 10 dB. Therefore, the bias caused by the PLL phase-tracking error 
will be less than 0.434 dB. 

A. 2. 5. 2 Amplitude-Estimate Variance 

The variance of the estimate a*(t) can be gotten as the sum of the variances 
caused separately by the noise and by the phase error. As can be seen from 
Fig. A— 7, the variance due to the noise is just the (constant) power spectral 
density of the noise times the noise bandwidth of the filter. The latter is 
given by 


OC 

B C " 27T I 

Jo 


|K C(j W )| du 


i. f a. 

' 21 JL i * - 


4x 


(57) 


AGC 


AGC 


and, as noted above, the former is given by 75.4 N Q /A . Therefore, the noise 
component of the variance of a*(t) is 


75.4 N q B c 


18.9 N 0 
A 2 t 

AGC 


= 75.4 


V_L 

L A 2 J 


(58) 


where in the last expression the term in brackets is the reciprocal SNR in the 
PLL loop bandwidth B , and is less than 0.1. 

Li 
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The second contribution to the variance of the estimate a*(t) is that due to 


the phase error ( t) . If the phase-error input of Fig. A-7 is again approximated 
by 

-20 log cos <}>(t) = 10 (log e) <£ 2 (t), (59) 

then to compute the component’s contribution to the variance, the power spectrum 

2 

of <f> (t) must be estimated, multiplied by the magnitude squared of the filter 

transfer function, and integrated over frequency. This can be done approximately 

o 

by assuming that the spectrum of <}>(t) is flat with a value of a /B from 0 to the 

L 

2 

PLL loop bandwidth B . Then the spectrum of <j> (t) is triangular with a peak 

L 

4 

value of 2o /B at zero frequency, decreasing to zero at a frequency of ±2B 
L L 

(see Ref. A5) . A good approximation to the variance at the output of the filter in 

Fig. A-7 can then be gotten by integrating this spectral density from 0 to the 

filter bandwidth B . This is equivalent to replacing the RC filter by a flat 

filter of the same noise bandwidth. The result for B < B is 

C _ L 



(60) 


Therefore, the variance in a*(t) due to the phase error <{>(t) is about 


2 4/V 

2(10 log e) a ( — 

\\j 


1 - y ■ 


Vl 


L A J 


2 /B_ 


1 - 


4B 


(61) 
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From the last expression and the one obtained above for the variance con- 


tribution due to noise, we obtain the desired result: 


var[a*(t)] = 75.4 


Vl 

. A 2 J 


1 + 


Vl 


L a 2 J 


l - 


4B 


L J 


(62) 


Since normally the SNR A /N B in the PLL loop bandwidth is at least 10, and the 

U ij 

bandwidth ratio 


B 4B x 
L L AGC 


(63) 


is smaller than unity, this variance is well approximated by 


var[a*(t)] ~ 75.4 


Vl 

L a 2 . 


(64) 

2 


which is just the component due to noise, and which will be less than one (dB) 
for x > 2/B . In other words, the contribution of the phase error to the 

AbL Ij 

variance is normally negligible, so the approximations we used in calculating it 
are not critical. 

A. 2. 6 Amplitude Tracking 

Since subsequent processing of the amplitude estimates a*(t) will be aimed 

at determining the power spectrum of the fluctuations in a*(t), the most critical 

information required about the amplitude tracking (estimation) errors in a*(t) is 

their spectrum. As shown above in Section A. 2. 4 and Fig. A-8, at frequencies well 

below 1/2 itt the spectral level of the noise-induced error in a*(t) will be 
Abb 

2 2 

75.4 NqA (dB) /Hz. At higher frequencies this will be modified by the squared 
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magnitude of the equivalent RC filter of Fig. A- 7, and will fall off by factor of 


100 for every decade increase in frequency above 1/2 ttt . 

Abu 

Another important consideration is the selection of the time constant t 

AGC 

Clearly, this should be small enough to pass the actual fluctuations of interest 
in the true signal amplitude a(t) . However, if it is too small the noise induced 
error in a*(t) will cause the amplitude tracking error a*(t) will cause the 
amplitude tracking error a*(t) - a \^-j l.g ue loo xarge iur tne linear model or 


Sections A. 2.4 and A. 2.5 to be valid. 

The approximation which leads to the linear model is replacement of the 
signal component 


— cos (0 1 - e 2 ) - i 
g 


(65) 


at the input of the integrator of Fig. A-5, by 


Jin 


[ir cos (6 i ' V] 


( 66 ) 


This is a valid approximation so long as 




Jin I cos ^ 6 i “ 0 2^ I <<2 
g 


■] 


(67) 


or 


20 log A +20 log cos (0, - 0 ) -20 log W <<17.4 

12 g 


( 68 ) 
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which is the same as 


a(t) - a*(t) +20 log cos (0 - 0 2 ) «17.4 dB 


(69) 


We have already seen in Section A. 2. 5. 2 that the standard deviation of the left 
side of this relation is well approximated by 

' 1/2 


jvar[a*(t)] W 2 = . 

18.9 N„ 

0 

2 

1 J 

A 

, AGC , 


(70) 


Therefore, for reliable tracking we must have 


18 - 9 V 

I a ^ t agc . 


1/2 


«17.4 . 


Taking "«" to mean less than 10%, this can be written as 


t agc 


18.9 N 


(1.74A) 



6.24 N q /A 



(71) 


(72) 


This condition was satisfied for all data contained in this report. 

A. 3 SPECTRUM ANALYSIS OF THE AMPLITUDE ESTIMATES 

Once the amplitude (in dB) as a function of time has been estimated as 
described in Section A. 2, it is possible to separate this estimate into three 
components: (1) the fluctuations due to Venusian atmospheric turbulence, (2) the 

(slowly varying) "trend" component due to average atmospheric attenuation, antenna 
beam-pattern effects, and changing earth-spacecraft distance, and (3) the varia- 
tions due to noise at the receiver input. The separation is made in the spectral 
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domain, and is done in two steps. First, the amplitude estimates a*(t) are 
"detrended" by fitting a quadratic time function to them, and subtracting it from 
them. Then, the difference is spectrum analyzed using an FFT algorithm, and 
theoretical results are used to identify the turbulence components and the noise 
components. 

A. 3.1 Amplitude Sampling 


14 0 -i-l 


As described in Sections A.l and A.z, tne processing which jj.cj.uo cue 
normalized amplitude estimates a*(t) is performed digitally. The time sampling 
rate of this digital processing is determined by the bandwidth of the processed 
signal (see Section A. 1.4). However, the amplitude estimates produced are low- 
pass filtered by the response of the (digital) amplitude estimating circuit, as 
shown in Fig. A- 7, to the spectral shape shown in Fig. A-8. Though these esti- 
mates are not strictly band limited, at frequencies above 1/2 t agc their spectrum 

2 

is attenuated by more than a factor of 1 /tt = 0.1. Therefore, spectral components 

above a frequency of 1/2t might be neglected, so that samples of the amplitude 

AGC 

could be taken at a rate as low as 1/x , which corresponds to a sampling period 

Auu 

as long as Spectral components above a frequency of 1/2 t agg above 


AGC 


2 ttt f = tt in Fig. A-8] would then be "folded" down below this frequency, but 

AGC 

only those above a frequency of 


2t 


AGC 


2t . 


AGC 


2ttt 


AGC 


= M1 


/2t 


AGC 


would be folded below a frequency of 1/2 itt , and these would be attenuated by 

Aul/ 


more than a factor of 1/ (2 tt- 1) = 0.036 =t -14dB. Thus, if the amplitude estimates 
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a*(t) are sampled at a rate no lower than 1/x , the spectrum of these estimates 

AGC 

will exist up to a frequency of l/2x , and all frequency components below 

AbL 

1/2ttt^ g c w iH be essentially undistorted. 

In the processing described in other parts of this report, the amplitude 

estimates were sampled at a rate no lower than l/ T APr ,« 

Abb 

A. 3. 2 Normalization and Detrending 

The theoretical turbulence results contained in this report are expressed in 
terms of the normalized amplitudes x(i)» which are defined as 


X(t) = In 


A(t) 

A(tj 


(73) 


where 

A(t) = received signal rms amplitude, 

A(t) = ensemble average value of A(t) . 

This can be expressed as 

x(t) = In A(t) - In A(t) 

= ^ [20 log A(t) - 20 log A(t) ] 

= 0.115 [a(t) - a(t) ] (74) 

where 

a(t) = received signal amplitude in dB, 
d(t) = A(t) expressed in dB . 
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The &(t) term is slowly varying, and in the processing is lumped in with the 
other "trend" components. Thus, our estimate of the fluctuations in x(t) -*- s 

x *(t) = 0.115 [a*(t) - Quadratic Fit] (75) 

where 

a*(t) = estimated received signal power in dB, 
and the "Quadratic Fit" is obtained by making a minimum least squares fit to a*(t) 
over the time interval to be analyzed. This essentially removes all slowly- 
varying components from x*(t). 

A. 3. 3 Spectrum Analysis 

The samples of x*(t) which result from the detrending process, are processed 
through a fast Fourier transform (FFT) routine to estimate the spectrum of x(t) • 
This is a conventional FFT routine augmented by two special features : a Hamming 
window and zero fill. The characteristics of these features are described in 
detail in Ref. A6. The purpose of the Hamming window is to reduce the spurious 
"sidelobe" response of strong spectral components (with a corresponding loss in 
frequency resolution) , and the purpose of using zero fill is simply to evaluate 
the Fourier transform of x*(t) at closer-spaced frequencies than would otherwise 
be obtained. 

A. 3. 3.1 Hamming Window 

If we let h denote the sampling period used to sample x*(t) > and we compute 
the power spectrum (magnitude squared of the Fourier transform) for N equispaced 
samples of x(t), then the time interval processed is T = Nh, and an FFT will yield 
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l+N/2 samples of the spectrum spaced in frequency by 1/T Hz from 0 thru l/2h Hz. 

A single (sine wave) input spectral component at frequency f Q will produce an FFT 

2 

spectrum output proportional to sine [(f - f Q )T] where 

sinc (x) = s JlLJZ . (76) 

TTX 


As a function of frequency f, this function has a main lobe of width 2/T centered 

at f , surrounded by "sidelobes" of width 1/T, the first of which is 12 dB below 
o 

the main lobe, with subsequent sidelobes falling off at 6 dB per octave. From 
this it can be seen that if the data to be spectrum analyzed contain wide fluc- 
tuations (greater than 12 dB) in spectral level, then the sidelobes from the 
strong components can mask the weaker components. To prevent this, a Hamming 
window is applied to the data. This is accomplished by multiplying the input to 
the FFT by the window function 

w (t) = 1 - k cos (27rt/T) 0 < t < T = Nh (77) 

where 

k = 23/27 for Hamming window. 

When this is done, a single spectral component at frequency f Q will have a spec- 
trum with a main lobe of width 4/T (thus, poorer frequency resolution) centered 

f t surrounded by sidelobes of width 1/T, all of which are at least 40 dB 
o 

below the main lobe. Thus, use of the Hamming window prevents sidelobes from 
strong spectral components from masking weaker components that are no more than 
four orders of magnitude weaker. The cost for this is a somewhat reduced 
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frequency resolution. For the data analyzed in this report, frequency resolution 
was not a significant problem, but a large dynamic range was,. so a Hamming window 
was used for all spectral estimates. 

A. 3. 3. 2 Zero Fill 

The Fourier transform of x*(t) exists at all frequencies, but is evaluated 

by FFT routines only at frequencies spaced by 1/T Hz. Thus, interpolation is 

required to get a smooth curve of the spectrum versus frequency. An efficient 

way to accomplish this is to use zero fill. In doing this, the original N data 

points (weighted by the Hamming window) are followed by N zeros, and the FFT is 

z 

applied to the resulting (N + N ) -point sequence. It can be shown that the only 

effect this has on the spectrum is that it is evaluated every 1/(N + N )h Hz, 

z 

i.e., the spacing between evaluation points is reduced by a factor of 1 + N /N, 

z 

producing a smoother representation of the curve. For the spectra of this report 
the values 


N = 1500, 

N = 2596, 
z 

N + N = 4096, 
z 

were used, which yield frequency points spaced by only 0.366/T Hz instead of 
1/T Hz. It should be emphasized that this procedure does not improve the frequency 
resolution , which is determined by the Hamming window, but only interpolates to 
provide a smoother curve. 
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A. 3. 3. 3 Spectrum Normalization 


For direct comparison with the theoretical results it is necessary to 
maintain the correct spectral level at the output of the spectrum-analysis 
processing steps. This means that the effect of the Hamming window must be 
accounted for, as well as that care must be used in determining the appropriate 
normalization at the other steps. 

The Hamming window is accounted for by computing the mean-squared value of 
the data x*(t) both before and after the Hamming window is applied. Then the 
spectral levels produced by the FFT routine are increased by the ratio of these 
mean-squared values. This assures that the area under the spectral curve equals 
the energy in the data, as is required by Parseval's theorem. 

A. 3. 3. 4 Frequency Averaging 

The final step in producing the experimental spectra for this report is 
averaging over frequency. This step is required due to the fact that we are 
estimating the power spectrum of a random process (the turbulence-induced ampli- 
tude fluctuations) from one sample function of that random process, and estimates 
obtained this way are very noisy. For example, the standard deviations of such 
estimates are approximately equal to the mean spectral level. This is illustrated 
further in Section A. 3.4. By averaging several adjacent frequency points together, 
the variances of the estimates are reduced, at the expense of reduced frequency 
resolution. 

Frequency resolution is most important at the lower frequencies, since 
frequency is most naturally shown on a logarithmic scale for results of the type 
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derived here. The lowest frequency which can be shown on such a scale is equal 
to about half the frequency resolution. Thus, the tradeoff to be made is between 
accuracy of the measured spectrum and definition of the low-frequency portion of 
the spectrum. 

The averaging done here is a "running-point" average, i.e., the spectrum at 
every frequency point (spaced by 1/ (N + N^)h Hz) is replaced by 



(f) 


1 

m 


(m-l)/2 

G(f + tAf) 

£=-(m-l)/2 


(78) 


where 

m = number of points (assumed odd) in running-point average, 

Af = frequency spacing from FFT, 

= 1/N + N )h, 
z 

G(f) = spectrum from FFT at frequency f. 


For the case m = 1, no averaging is done, and the frequency resolution is 
just that of the Hamming window. This is given by 


B 

w 


/» 00 


V —00 


F{w(t) } 



F{w(t) } 


f=0 


(79) 
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where F denotes the Fourier-transform operation, and 
w(t) = Hamming-window weighting function, 

= 1 - k cos (2?rt/T) , 

T = time duration of data spectrum analyzed, 
= Nh. 

The Fourier transform of w(t) is 


F{w(t)} = — e -^ 7tT ^ sin (irTf) ~ + 


kf 

2 2 
(I/O - f Z J 


so that 


L (I/O - f J 


and 




sin (irTf) 
7T f 


T. 


(80) 


(81) 


(82) 


This is the square root of the denominator of the expression for B . The 

w 

numerator is most conveniently evaluated by using Parseval's theorem, viz.. 



| F{w(t) }| 2 


df 




dt 


[1 - k cos (2Trt/T)] dt 


= T (1 + k /2) . 
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Therefore, we have for the frequency resolution (i.e., for the bandwidth) of the 
Hamming window, 


B 

w 


T(1 + k 2 /2) 

2 

T 


(1 + k 2 /2)/T 


| [1 + (23/27) 2 /2]/T = 1.36 3/T for the Hamming window 

|l/T for no window (k=0) . (84) 


As previously indicated, adjacent frequency points from the FFT are 
separated by 


Af 


1 

(N + N )h 
z 


1 

(1 + N /N)T 
z 


1 0.366 

(1 + 2596/1500)T T 


(85) 


Therefore, when m points are averaged together the frequency resolution becomes 

B = B + (m - 1) Af 
w 

= [1.363+0.366 (m - 1)]/T. (86) 

The lowest frequency point of interest (i.e., above zero frequency) after this 
averaging is 
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f = Af = 0.366 (m + 1)/2T 

min 2 

_ B 0.316 
2 T 


( 87 ) 


A. 3. 4 Confidence Limits 

It is shown in Ref. A6 that the spectrum estimates we obtain are approxi- 
mately chi-square distributed with a number n of degrees of freedom given by 
(Ref. A6, eq. 6.110) 

n = 2BT ‘ (88) 


where 

B = frequency resolution, 

T = time duration of data spectrum analyzed. 

For our parameter values this gives 

n = 2.726 + 0.732 (m - 1) . 


It then follows that a 90% confidence interval is given by 


n 


n;0.05y 


g av <£) £ 


n;0.95 , 


g av <£) 


(89) 


JPL Technical Memorandum 33-644 


A- 37 



where 


G (f) = spectrum estimated from running-point average, 

G(f) = true spectrum, 
o 

y = percentage points of the chi-square distribution. 

n;a 

Values of y 2 for integer n are tabulated in Table A. 3 of Ref. A6. These were 
n ; a 

used to tabulate the 90% confidence intervals shown here in Table A-l. 

From that table it can be seen that without averaging, the 90% confidence 
interval is 

0.38 G AV (f) < G(f) <8.5 G AV (f) 

which is a very wide interval. When a 51-point average is used this interval is 
reduced to 


0.72 G AV (f) * G(f) < 1.5 G AV (f) 

which corresponds to a much more useful estimation accuracy. 

A. 4 SUMMARY OF PROCESSING STEPS 

In this Section we will summarize the processing steps required to produce 
the experimental spectra shown in this report. 

A. 4.1 Filtering and Decimation 

The data processed for this report had previously been received, recorded, 
and digitized for the original Mariner V occultation experiment (Ref. A7) . Thus, 
the first step to be taken here is to filter and decimate the data (Section A. 1.4). 
This requires selection of a spline curve to follow the gross frequency-versus-time 
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profile of the received signal, and a choice of bandwidth A and decimation factor 


N. The spline curve had been previously developed, as described in Ref. A8. The 
bandwidth A must be chosen large enough to pass the signal frequency variations 
after being shifted by the spline curve, and the decimation factor N must cor- 
respond to A. For the processing carried out here, N = 101 was found to be a 
convenient value. 

A. 4. 2 PLL Tracking 

The main parameters to be chosen for the phase-locked- loop tracking (Section 
A. 2) are the loop bandwidth B , the damping R, the AGC time constant x , and 

Li AuL 

the output sampling period h. The bandwidth must satisfy 


4r 


< B. 


10 N, 


where 

R = signal's frequency rate (Hz/sec), 

A = signal rms amplitude, 

Ng = noise one-sided spectral density (watts/Hz). 


(90) 


The choice of damping R is not critical, and here was taken to be 2 (see Fig. A-4) . 
The time constant x is one of the more critical processing constants 

Abv 

(Sections A. 2. 5 and A. 2. 6), because it determines the maximum frequency at which 
amplitude fluctuations may be observed. For good performance it must satisfy 

t agc V a2 - <91) 
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Also, it must satisfy 


_1 < - (92) 

2tt f T AGC 

max 

where 

f = maximum frequency of turbulence, 
max 

One other consideration in the choice of x is that it indirectly determines 

the lowest frequency component that can be measured. This is because the lowest 
frequency component measured is (Section A. 3. 3. 4) 


_ 0.366 (m + 1) 
min 2T 

where 

m = number of points in running-point average, 
T = Nh, 

N = number of data points in input to FFT , 

= 1500, 

h = sampling period, 


(93) 


Thus, for fixed values of the other parameters, the minimum frequency f . is 
’ mm 

limited by the choice of t.__. 

AGC 

A. 4. 3 Spectrum Analysis 

The spectrum analysis step consists of fitting a quadratic function of time 
to the data for use in detrending it, computing the FFT of the Hamming-windowed 
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data, and smoothing the result by the running-point average. The only parameter 
to be chosen here is the number m of points over which to compute the average. 

In practice several values of m are used, usually ranging from one to 51. 

A. 4. 4 Interpretation of Results 

The last step in the processing consists of using the theoretical results of 

this report to separate the turbulence-induced fluctuations from the noise-induced 

fluctuations. The noise is recognized (see Section A. 2. 5. 2) by having the spectral 

2 2 

shape shown in Fig. A-8 and a level given by 75.4 N^/A for a*(t), and by N^/A 
for x*(t) , where a*(t) is the estimated amplitude in dB, and x*(t) is the natural 
logarithm of the estimated amplitude. The part of the experimental spectrum 
which deviates from this noise part is then interpreted in terms of turbulence- 
induced phenomena. 
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No. of Freq. 
Points 
Averaged 


n = 2BT 

Nearest Integer 
in Table A. 3 
of Ref. A6 

n 

Degrees 

of 

Freedom 

2 

X n;0.05 

2.73 

3 

0.38 

10.05 

10 

0.55 

17.38 

17 

0.62 

24.70 

25 

0.66 

32.03 

30 

0.68 

39.35 

40 

0.72 


n;0.95 


Table A-l. 90% Confidence Intervals for Averaged Spectrum Estimates 
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Fig. A-5. Block Diagram of the Coherent Amplitude Estimating and AGC Computer Program 
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Fig. A-6 . Block Diagram of the AGC Loop Analyzed by Tausworthe. Taken from Fig. 7-1 of Ref. A1 
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APPENDIX B 


EFFECTS OF ATMOSPHERIC TURBULENCE ON THE RADIO SIGNAL OF A 

FLY-BY SPACE PROBE 


B.l Formulation of the General Problem 

Let us consider both plane and spherical wave propagation with the respective 
configurations depicted in Figs. B-l and B-2. The direction of wave propagation is 
assumed to be along the x axis and the turbulence medium characterized by the 
refractive index 


n(r, t) = 1 + n^ (r, t) (1) 

where n^(r, t) is the fluctuation and assumed small, r is position and t time. 

According to Rytov's method (Refs. B1 and B2) the field U at the observation 
point (L, p ) at time t is given approximately by 

U(L, p, t) = U (L, p) exp [* (L, p, t) ] (2) 

o I 


^(L, P, t) - 




x* , p’ - yp) n^x, p' , t) 


(3) 


where 


h(L, x', p') 


2 r -*■ 

k 1 . k p' • p 1 

2ir y(L - x’) 6XP L 1 2 y(L - x')_ 


y = 1 for a plane wave 


Y 


x 

L 


for a spherical wave 
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p =yy+zz, y and z are unit vectors in the y and z directions, 

_ y» y + z' z, k is the free space wavenumber and is the field in free 

space (without turbulence). Also, 

^(L, P, t) = X(L, P, t) + i S 1 (L, p, t) (4) 

where X is the log-amplitude fluctuation and is the phase fluctuation. 

Under Taylor's frozen-in hypothesis 

n^r, t) = n 1 (r - vt, 0) (5) 


where v is the average wind velocity perpendicular to the direction of wave 
propagation. Inherent in (5) is the assumption that the effect of wind velocity 
fluctuations is small compared to the average wind velocity v. Applying Taylor's 
frozen-in hypothesis to (3) , we have 


^(L, P, t) = 



- (yp - vt)] n^(x, p" , 0) (6) 


where the coordinate p" = p' - vt is attached to the moving turbulence, in 

(6) is a function of position only and the usual time-independent spectral 
representation for n^ may be used. 

We will assume that the turbulence is described by a random process studied by 
Silverman (Ref. B3) and employed in previous turbulence studies (Refs. B4 and B5) . 
The correlation function for the refractive index fluctuation B^(r^, r 2^ ’*" S ® : *' ven 
as a product of the average coordinate and a function of the difference coordinate, 

i.e. 


B ( V V = B n^ c } B n ( V 

n 1 2 n c n a 


(7) 
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where r = ~ (r + r„) and r = r, - r„. The corresponding spectrum $ (K. , K ) 

c 2 1 2 d 1 2 nl2 

is then a product of a function of the average wavenumber and a function of the 
difference wavenumber, i.e. 


$ (K 
n 


1 * 


KJ = * V (K,) $°(K ) 
2 n d n c 


( 8 ) 


where K = — (K + K ) is the average wavenumber and K, ■ 

c 2 1 2 6 d 

VO VO 

difference wavenumber. B and B are related to $ and $ 

n n n n 

Fourier transforms 


is the 

through the following 


B V (r ) 
n c 


exp(i K * r ) 4> (K ) dK 
d end d 


(9) 



r ,) $ (K ) dK 
d n c c 


( 10 ) 


Since the wave is propagating in the x direction we need use the spectral 
representation in the y-z plane only. We then have 



exp (i k 


P ) 
c 


F n ( V 


-V -v 

O dK. 


( 11 ) 


Vv p d ) ' 



exp(i k • p ,) F (x,. k ) tk 
c d n d c c 


—00 


( 12 ) 
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and the two dimensional spectra F V and F° are related to the three dimensional 

n n 

v o 

spectra $ and 4> by 
r n n 


F (x , k ) = 
n c d 


ex P <i MV V K dl" K d ) “dl 


(13) 


?°(X , K ) = / 

n d c I 
J _ 


exp (i K x ) $ (K , k ) dK 

cl d n cl c cl 


(14) 


where K = K x + K y + K z = K x + k , K = K x + K y + K z 
„ c ^ cl c2 c3 _^cl ^ c d dl d2 d3 

■ K dl X + V ' V + V 72, ^d " - V X c ‘ (X 1 + X 2 ,/i a " d x d ‘ X 1 - V 

We will assume that the turbulence is isotropic so that the correlation function 

B°(r ) is a function of the magnitude of r, only and (Ref. B4) 
n d d 


F°(x , k ) dx = $°(0, k ) = $°(k ) (15) 

ndc d n c nc 

The correlation functions for the log-amplitude and phase fluctuations, 

B and B , are 
X s’ 

B (L, p , p , t , t) = <x(L, p , t ) X(L, p , t )> (16) 

Xcdc 11 22 

B ( L » P , P , t , x) = <S (L, p , t ) S (L, p , t )> (17) 

s cdc 1 111 22 

where t^ = (t^ + t^)/2 and T = t jL "" t 2* Following a procedure similar to that 
in Ref. B4, we can obtain (16) and (17) from (4) and (6) for the turbulent 
medium described above. Let us first look at the case of log-amplitude fluctu- 
ations. We will assume that the observation distance is large compared to x 

d 

and obtain the following results: 



B-4 
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B (L, p , p , t , x) = 2nk 
X c d c 


r L x f[ - v - 

0 CO 


G (x , yp - vt - p , yp - vt) 
X c c c c d 


( 18 ) 


where 


If- 


G (x , p , p_) = I I dK exp (i k • p ) 
X c c d I f d d c 


// 


die exp(i k • p.) S (x , ic , K ) $ (k ) (19) 

c c d x c c d n c 


S (x , k , k ) = — {exp(i A k • K,) + exp(-i A k • k,) 
X c c d 4 I cd cd 


- exp 


2 d 
i A 1 k + ~T 
c 4 


- exp 


2 


-iAU 2 + - d 
c 4 


f (20) 


and 


A = 


Y(L - x ) 
c 
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II -*■ 

We note that B V (x , p ) is measured in the coordinate system p £ which is moving 

with velocity 5, i.e.! ?' = p' - vt where p’ is the actual position of the tur- 

bulence at t , We can transform (18) into an integral with respect to the actual 
c 

f a 

coordinate P and obtain: 
c 







27Tk 




• G (x , YP " p', YP. " vr) (21) 

X c c c a 

It should be pointed out that (21) can also be written as an integral in terms of 

the two wavenumbers tc and k as was done in Ref. B4. 

c d 

Similarly, for phase, we have 


OO 



—00 


G_(x , ^ p - P > " VT > 

See c d 


( 22 ) 


where is the same as (19) except for S which is replaced by S and 


s s ( v V K c> ■ I 



/ 2 K d\ 

exp(iA< • k.) + exp(-iAK • k ) + exp 
c d c d 

_ iA \ K c + ry 


+ exp 


2 v n 

2 <d 
-iA l k + ” 
c 4 


(23) 


B-6 


JPL Technical Memorandum 33-644 



2 J 2 

The variances of the log-amplitude and phase fluctuations a and a 

X 15 

respectively along with the frequency spectra of the log amplitude and phase 
fluctuations W and W^ respectively may now be derived in terms of the corre- 
sponding correlation functions (18) and (22) . Since we are interested in the 

signal at one location and at one time, p, = 0 and x = 0. 

d 


a (L,p ,t ) = B (L, p , 0, t , 0) 
X c c x c c 


(24) 


a Q ( L > p .t ) = M L . p , o, t , o) 

S c c S c c 


(25) 


and 


/ 


W (L,p ,0,t ,f) = W (L, p , t , f) = 2 I B (L, p , 0, t , x)e lt0T dT (26) 
xcc Xcc fXc c 


/ 


W c (L,p , 0,t ,f) = W (L, p , t , f) = 2 I B (L, p , 0, t , x)e ^dx (27) 
See bcc |Sc c 


where co = 2irf and f is the fluctuation frequency. We will evaluate and discuss 
(24) -(27) in the following sections. 

B.2 Theoretical Variances of Log-amplitude and Phase Fluctuations 

As in a previous study (Ref. B5) we will assume that the Venusian atmospheric 
turbulence is localized and smoothly varying so that its correlation function 
can be written 


B 

n 




0.033 c 2 (r ) 
n c 



(28) 
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where the structure constant c characterizes the intensity of turbulence and 

n 2 -*■ 

is slowly varying. Comparing (28) with (7) we observe that B n ( r c ) = °* 033 c n ^ r c ^ * 

We will also employ the von Kdrmdn representation of the Kolmogorov spectrum to 

describe the localized homogeneous and isotropic turbulence, i.e. 



(29) 


which is of course related to B (|r |) through (10), (14) and (15). 

n d 

Using this turbulence description, G in (18) can be integrated so that (24) 

X 

and (25) become (Ref. B5) 


2 -*■ 
a (L, p 
X 


t ) 
c 


= 0.033 



00 




-V.II 

p 

c 


- vt ) 
c 



/ 
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where Re [ ] means real part of [ ] and 4>(a,c;z) is the confluent hypergeometric 
function (Ref. B6) . 

As discussed in the main text (Section 4.2) we will choose the following 
2 

model for c 

n 


2 2 2 
c (r ) = c (x ,y ,z ) = c exp - 
n c n c c c no 


M ■ ft? - ft)l 


(32) 


2 2 

where c = c (0,0,0) is the peak structure constant and the coordinate system 
no n 

and occultation geometry are shown in Fig. B-3. 

It is obvious that the moving spacecraft results in a sweeping motion of the 
line-of-sight through the turbulent atmosphere and that this motion has the same 
effect on the radio signal as the wind. Also, since the velocity of this sweep- 
ing motion is considerably larger than the wind velocity in the Venusian 
atmosphere, v essentially represents the velocity of the sweeping motion and is 

oriented in the y -direction. It can therefore be seen that this study does take 
c 

into account the effect of the spacecraft velocity. 

Applying (30) and (31) to the configuration in Fig. B-3 we have at the obser- 
vation point located on earth 


= 0.033 — 

2tt 


f OO 

1 “II 


dy dz c (x ,y - vt , z ) 
c c n c c c c 


G (x ,y ,z ) 
x c c c 

G q ( x „>y„» z ) 

O c c c 


(33) 
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where 



( 34 ) 


(35) 


A 


x (L 1 - x ) 
c c_ 

kL’ 


and v the velocity of the I Sweeping motion of the line-of-sight is assumed to be in 
the direction. We note that A can be alternately written 


It 

x L 
. c 

ti 


k(L + L ) 


and since L » L, 


A can be approximated by x / k. 

c 


Furthermore 




B- 1 0 
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and when L >> a 



f (x ) exp 
c 



n/tt a f(L) 


(37) 


It can then be shown that (33) becomes 


2 

a 

X 


0.033 — n/tt a c (I + I ) 
2ir no 1 2 


(38) 


2 

a 

S 


k 2 

0.033 — nTtt a c (I - I ) 
2tt no 1 2 


(39) 


where 



(40) 
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As shown in Ref. B5, (38) -(41) can be evaluated for various ranges of L. 
In this report we will be interested only in the first three regions which are 
shown in Fig. B-4. 


L 2 2 2 

Region I; A = << L << L b << L a << b << a 

° k o o o 


We will assume that A/L << vt , so that 

o c 




2 4 . 2 

y + z 

i i c c i 

d2 ^ A 2 V 

—CO ' o 


, -11/6 

I 

< 

rt 

O 

to 

1 

1 exp 

L‘ J 


(42) 


12 3 5/3 

= — tt L exp 
5 o 




(43) 


Also, using the small argument approximation for \p (Ref. B6) 




5/6 


(44) 


where T is the Gamma function, we have 


*2 = 


, 2 5/3 

2f L TT 
o 


6 , „ 

( 5\ 

/5tt\ 

— + r 

- — 

cos 1— r) 

5 

\ 6/ 

\12/ 

_ 




£11 

5/3 


o J 


exp 


-(=*)“ 1 


(45) 


so that 


2 , 7/6 5/6 2 

a =0.99 7k L ac exp 
X no 


-m 


(46) 


B— 1 2 
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2 5/3 2 2 \ 

,385 k L a c - a ) exp 
o no x / 


where a in (47) is just (46) 
X 


2 2 2 
Region II: L « A « L b « L a « b « a 

o o o 


In this region I is the same as (46) while for I we use the large argument 


approximation for ^ (Ref. B6) 




so that I^ becomes I^ = 0 and 


2 2 „„„ k r 2 12 3 5/3 

a = a = 0.033 — \rrr a c ~ n L exp 
X S 2tt no 5 o 


2 5/3 2 

= 0.69274 k L a c exp 
o no 


; w; 




2 2 2 
Region III: L << L b << a << L a << b << a 

o o o 


. / 2 \ -11/6 r 2 2 " 
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and can be rewritten 



( 51 ) 


(52) 


I is the same as that in region II. Therefore 


bL 5 ^ 3 r 

a 2 = a 2 = 0.033 c 2 a ir 2 k. 3 — r 

y S no L 


fl) 


ft) 


exp 




(53) 


T 8/3 

3 2 

0.548 k ab ° 


~ c exp 
L no 


h t \2 

c 


(54) 


It should be emphasized that we have used the approximation A/ << vt^ in the 
above results and that this condition should be checked especially in region III. 
We note that the results contain the factor 


exp 



which becomes significant only when vt^ > b. This factor represents the effect of 

the finite extent of the turbulence in the y direction when the line-of-sight 

c 

traverses a distance longer than b. 


B-14 
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2 

The dependence of o on L is summarized in Fig. B-4. If vt^ << b , then 

2 X _ 

the behavior of o can be explained in a similar manner as in Ref. B5. In region 
X 

I we see only the effects of turbulence inhomogeneity and finite extent in the 

x direction, while in regions II and III the effects in the transverse direction 
c 

as well. In region II the largest blob size is L so that the scattering angle is 

—1 2 ° 

no smaller than (kL ) . Also since L » kL corresponds to the Fraunhofer region 

° ° -2 

of L , the contribution from each blob decreases as L . Therefore, even though 
o 

the n umb er of blobs contributing to the fluctuation increases as the square of 
the distance L, the variance of the log-amplitude fluctuation remains essentially 
constant. 

At L = kL b, the Fresnel-zone size is the same as the size of the turbulent 
o 

medium in the y direction of the L = 0 plane. In region III, L > kL b; the 
c ° 

number of blobs increases as L and since the contribution from each blob decreases 

—2 # . -1 
as L , the variance decreases as L 

2 2 

It should be point out that even though the expressions for a ^ and o g in 

regions II and III are asymptotic, they are also valid close to the boundary 

separating regions II and III. We can obtain some estimate on the error in using 

2 

these results close to the boundary by evaluating and comparing a at the 

X 

boundary L = kL b using (49) and (54) . We find that the ratio of (49) to (54) at 

L = kL b is 1.26. 
o 

B. 3 Theoretical Frequency Spectra of Log-amplitude and Phase Fluctuations 
Let us rewrite (26) in the following form 


00 



—00 


(55) 
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where 







00 



_oo 


exp(iK • vt)S 


(x ,K.,K ) 
c d c I 



n ' c' 


-itox , 
e dx 


(56) 


and S is given by (20) . 
X 


For the x integration in (56) we use the fact that 


1 _ 
2t r 



vx - itox) dx = S (to 


-y 
+ K 


V) 


(57) 


where 6 is Dirac's delta function so that 


®X ( 


x 


-J f d * d (w c -^ j 


■y , -y ■>. . -y . o , . 

d< 6 (to + k • v)S v (x , K j k )$ (k ) 

c c Xcdcnc 


(58) 


For the k integration we use a polar coordinate system where k is the radial 
c c 

distance and a is the angle k makes with v. We will first consider the integra- 

c 

tion over a. Since (Ref. B7) 



da 6 [f (a) ] = 



(59) 
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where a are the roots of f(a) = 0 within (a, b) , we obtain 
n 


I. 


[ ]/( K.V ) 2 - J 


daS(u> + k v cos a) = < 

c 1 


d a [6(a “ c^) + «S(ct- o^) ] for k^v > u) 


(60) 


for k v < a) 
c 


where 


a 


1 


-1 

cos 



a 


2 


a 


1 


0 - - it 


We denote the vectors k corresponding to a, and a., 

c 12 

(see Fig. B-5) . Then (46) becomes 


k and k „ respectively 
cl c2 






exp[iic , 
d 




—CO 



1 

K dK 



[S (X ,K ,K ) + S (X ,K ,K )] $ (K ) 
X c d cl x c d c2 n c 


(61) 


Let us now carry out the integration over 

d 

that the wind direction is along the y axis. 

c 


Before we do this we will assume 
We will also write 
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B = YP„ “ P 
c c 


C 1 ■ 


c 2 - 


Note from Fig. B-6 that C, = C„ = - A co/v and C, = -C , where the subscript 

ly 2y lz 2z 

c c c c 


y or z designates the 
c c 


icopcCtlVc Cuuijjuucu L • O J-Iice 


00 

JJ • 


B + iA< + i 
c 


. A 2 \ - 1 - 4 it / 2 tt B^ | 

L 4 K d/ dK d = r exp ( iAK c + i 2' i Tl (63) 


c 2 A 


(61) becomes 


; (x , p , p ) = 
X\ c c c / 


k dK 

c c 

I 2 2 

'(k v) - fa) 

c 


*v ) 

n c 


IM 2 


f (B + C ) 6(B + C, ) 

I y ly z lz 


+ {<B y - C ly > S(B Z - C lz > + 6 « y + C ly > «<B - C ) 

c c c c y c Ay c C 1 C 

+ 6(B y c ‘ C ly c } 5 (B Z c + C lz c ) ]' A 1 Re 6Xp ( iA< c " 1 ^2 + 1 2 ) (64) 

where Re [ ] means real part of [ ] . For the integration over k , we will use 

c 

the following change of variable 


2 w '2 

K - = K 

c . 2 c 

v 
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and from Fig. B-6 we note that C 


= -C = Ax' = Ak ' . Thus, 
lz 2z cl c 

c c 


g ) 

X c c c 


dK 


$°(k') f 6 (b — A - Wb + Ak )+ s(b + A -) 6(B - Ak ) 

n c 4 l \ y v z c \ y v/ z c 


*'0 


+ 6(B -A-J6(B - Ak ) + 6(B + A “ I 6 (B + Ak ) 

\ y v / z c \y^ v/ z c 


»|B + A - ) 6 (B + Ak )l 

\ y v/ z c i 


4ir 

— Re 
A 


■ . . 2 . B_ . tl 

exp | iAk c - x a + i 2 


( 66 ) 


and finally 


g x ( V 0 ■ S X1 ( V V *c> + *X2 ( V V 


(67) 


where 


®xi^ x c’ V ^ - ^2" | 6 ( B y - v) + 5 (v A v) 


p 


( 68 ) 


^2 






+ i A 


2 

U) 

2 

v 



(69) 


The frequency spectrum of the log-amplitude fluctuations for a turbulent medium 
characterized by (7), (8) and (15) is given by (55) along with (67). 


JPL. Technical Memorandum 33-644 


B- 19 



Following a procedure similar to that described above the frequency 

spectrum of the phase fluctuations W (L, p , t , f) can be derived. 

l) c c 


W (L , p ,t ,f) 
S c c 


o 2, 2 

8tt k 


r r 


dx 


v -*■' 
dp B (x ,p 
c n c c 


vt ) g (x ,p ,p ) 
c S c c c 


( 70 ) 


where 


8 S ( V V ■ 8 xl ( V V - g X2 ( V V 0 


( 71 ) 


B.4 Theoretical Frequency Spectra for Mariner 5 

V o 

The results given in Appendix B-3 were in terms of B and $ . In this sec- 

n n 

tion we will apply them to the case of Mariner 5. As in Appendix B2 we will 

assume B V and $° to be given by (28) and (29) respectively. We will also 
n n 

assume that the observation point is located on the x axis (Fig. B-3) so that 

-v ' c 

p = o. Then according to. (62) B = - z and 
c z c 

c 



o 


Substituting (72) into (69), 


g „ becomes 
b X2 


(72) 


S X2 


(x 

c 




exp 



2 

u 


2 

v 
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Let us apply (55) to the configuration in Fig. B-3. At the observation point 
on earth 


where 


W (l\ t , f) = 8 A 2 I dx [M (t , f) - M (t , f)] 
X c | c 1 c 2 c 

0 


(73) 


M (t f) = 
1 c 



dp B (x , p - vt ) g (x , p ) 
c n c c c XI c c 


(74) 


M„(t , f) = 
/ c 


dp B (x , p 
c n c c 


vt ) g Y (x , p ) 
c X2 c c 


(75) 


Then, using (68) and (73) and model for c in (32), we obtain 

n 


M, (t ,f) = 0.033 c exp 
1 c no 


(X - L) 2 
c 

v? /si A" 

1 1 A 2 /J l\ 

2 

a 

*V V 

2 '" s '. 2 ^v 2 L 2 J 



r /. a) , \2 1 


r /. co \ 2 i 








\ v c/ 


V v c/ 


exp 

, 2 

+ exp 

2 



b 


b 



(76) 
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M (t ,f) = 0.033 c 2 
2 c no 



To carry out the integration over x in (73) we note from Fig. B-3 that 

c i( 

A = (L - x ) x /kL can be written alternately as [x L / (L + L )]/k. L is 
c c c 

the distance from the planetary atmosphere to the observation point on earth 

and is much longer than L; A can therefore be approximated by x /k. Also, since 

c 

L » a, we can extend the lower limit of integration (73) from o to - 00 and 
noting that 
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and 


f 

(x - L) 2 x 2 

C . . C 

I dx exp 

2 + k 2 

L c 

a v 


f(x ) 
c 


= n/tt a exp 


i 


L_ uj_ 
k 2 


v 


1 

4 



f(L) 


if f(x ) is slowly varying, (73), (76) and (77) become 


(79) 



JPL Technical Memorandum 33-644 


B-23 



Similarly, 


W (l', t , f) 
S c 


0.033 c 2 8iT 2 k 2 

no 2v 



( 83 ) 


We can check if (80) and (83) reduce to the results for a plane wave 
propagating in a homogeneous turbulent medium by letting a -* 00 and b -* “ and 
noting that (Ref. B6) 


lim Ha, c; z) = ” ^ c ) if Re ^ < ° 

z+o 

where F is the gamma function. We obtain at = L 

W (f) = 8 irk 2 L (0.0 33) c 2 - [W (f) + W (f) ] (85) 

X no v XI X2 

W (f) = 8 irk 2 L (0.0 33) c 2 - [W v1 (f) - W (f)] (86) 

S no v Xl X2 


where 


V» 


\fn 

2 


(t) 




2 2 
L 

o , 


-4/3 


(87) 


W X2 (f) 


n/tt 


Re 




1 

2 ’ 


1 

3 ’ 


- l 


L / co_ 
k ( 2 


v L 


V 

2 

V J 


( 88 ) 


which agree with those derived by Tatarski (Ref. Bl) . 
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Let us evaluate (80) - (83) for regions II and III as defined in Appendix B2. 

2 

We recall that in regions II and III we have the condition that L q << L/k « L^a. 
We can therefore use the small argument approximation for ip in M* and the large 

i 

argument approximation for <J> in M 2 , i.e. (Ref. B6) 




We also note that since b » L , L/k « b in regions II and III and we finally 

o 

obtain 


W (f) = K 
X 




(91) 



(92) 


where 


8/3 


K = 0.033 c 2 4ir 3 
no 


k 2 a 


( 93 ) 
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The validity of Taylor's frozen- in hypothesis for wave propagation in 
localized smoothly varying turbulence deserves some discussion. We will follow 
similar arguments used by Tatarski (Ref. Bl) in his discussion of wave 
propagation in homogeneous turbulence. 

In the case of localized smoothly varying turbulence, the fluctuations are 


mainly attributable to inhomogeneities in the medium whose scale size l is of 

the order L moving with velocity v. (Of course in the case of Mariner 5 v is 
o 

the sweeping motion of the line-of-sight.) The inhomogeneities exhibit intrinsic 


random velocity in addition to the constant wind velocity v and this random 

1/3 1/3 

velocity is of the order (el) or (eL^) where e is the rate of dissipation 

of the turbulent energy. The lifetime x of these velocities (time during 

v 


which these velocities remain constant) depends only on e and l and is of the 
1/3 2/3 

order e l . On the other hand, the time required for an inhomogeneity moving 


with velocity v to cover a distance l = 1 is of the order T - L /v. The 

o o 

frozen-in hypothesis will therefore be valid when x » T or 


v >> (eL ) 
o 


1/3 


(96) 
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According to Refs. B8 and B9, for the atmosphere of Venus e is of the order of 

10 cm sec . As shown in the main text L is of the order of 100 m. For 

o 

these values, v must be greater than 0.0215 m/sec. This condition is certainly 
met in the case of Mariner 5. 
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Fig. B-2. Spherical Wave Configuration 
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Fig. B-3. Mariner 5 Occultation Configuration 
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Fig. B- 


6. Vector Diagram for and 
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